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PURPOSE 


The  purpose  of  this  program  is  to  investigate  any  methods,  tech¬ 
niques,  or  equipment  uhlch  can  be  used  to  determine  atmospheric  wind  fields 
at  locations  remote  from  the  sensing  and  measuring  apparatus.  The  work  is 
to  be  perfonned  in  three  phases: 

I.  A  survey  of  the  field  to  ascertain  the  state-of-the-art  through 
a  literature  search  and  a  study  of  current  efforts  at  governmental  and  non¬ 
governmental  agencies  and  Institutions; 

II.  A  critical  analysis  of  work  done  in  the  field,  as  found  in 
Phase  I,  to  determine  the  capabilities  of  the  methods  for  remote  measurements; 
and 

III.  Ihe  design  of  experiments  to  verify  the  results  of  the  analysis 
in  Phase  II. 


ABSTRACT 


Analyses  were  made  to  determine  the  feasibility  of  three  proposed 
methods  of  remote  wind  measurement  using  (l)  scattering  from  natural  atmos¬ 
pheric  turbulence]  (2)  electromagnetic  scattering  from  acoustic  waves;  and 
(3)  infrared  tracking  of  an  artificially  heated  volume  of  air,  or  "bubble." 

Use  of  natural  turbulence  as  a  sensor  will  require  (l)  additional 
data  on  distribution  and  characteristics  of  turbulence  from  ground  level  to 
one  mile  altitude,  (2)  correlation  of  turbulence  motion  emd  wind,  and 
(3)  radar  state-of-the-art  improvement  to  provide  consistent  detection  and 
measurement . 

Remote  wind  measurements  by  microwave  reflection  from  acoustic 
waves  have  been  demonstrated,  but  additional  experimental  data  are  needed  to 
determine  maximum  usable  range  and  the  effects  of  turbulence  on  the  acoustic 
waves. 

Remote  generation  of  a  heated  bubble  of  adr  does  not  appear  feasible. 


-  vii 


I. 


PUBUCATIQNS.  LECTURES,  REPORTS.  AMD  CONFEREWCES 


Publications ;  None 
Lectures :  None 

Reports:  Monthly  Letter  Reports  Nos.  5,  6,  and  7. 

Conferences i 

On  6  November  1961,  a  conference  was  held  at  USASRDL,  Evans  Area, 
Belmar,  New  Jersey, with  the  following  persons  present:  for  the  Signal  Corps  — 
K.  C.  Steelman,  Contracting  Officer's  Technical  Representative,  SMI;  H.  H. 

Grote,  Project  Engineer,  SMI;  D.  A.  Delainger,  Chief,  Meteorological  Division, 
SM;  R.  M.  Marchgraber,  Physicist,  SMA;  and  Dr.  H.  J.  aufm  Kampe,  Chief, 
Atmospheric  Physics  Branch,  SMA;  for  Midwest  Research  Institute  --  H.  L.  Stout, 
Assistant  Director,  Engineering  Division;  P.  C.  Constint,  Jr.,  Head,  Electronics 
Engineering  Section;  R.  W.  Fetter,  Senior  Engineer;  and  Dr.  P.  L.  Smith,  Jr., 
Senior  Engineer. 

This  conference  was  held  to  discuss  the  results  of  the  first  phase 
as  presented  in  Report  No.  1,  and  the  program  for  Phase  II. Analysis.  After 
discussing  the  nine  possible  methods  suggested  for  analysis  in  Report  No.  1, 
it  was  agreed  that  only  three  of  these  methods  would  be  analyzed  to  determine 
their  capabilities  for  remote  wind  field  measurements.  These  methods  are: 

(l)  scattering  from  natural  atmospheric  turbulence,  (2)  electromagnetic  scatter¬ 
ing  from  acoustic  waves,  and  (3)  infrared  tracld.ng  of  an  artificially  heated 
"bubble"  of  air. 

On  19  and  20  February  1962,  a  conference  was  held  at  Midwest  Research 
Institute  with  Mr.  H.  H.  Grote,  Project  Engineer,  USASRDL,  and  Messrs.  P.  C. 
Constauit,  Jr.,  B.  L.  Jones,  auid  R.  W.  Fet+er  of  Midwest  Research  Institute. 

Principal  subjects  discussed  vjre  the  results  of  the  Phase  II  investi¬ 
gation,  the  preparation  of  this  phase  report,  and  the  evaluation  of  the  three 
methods  of  wind  measurement  analyzed:  .1)  electromagnetic  scattering  frexn 
natural  turbulence,  (2)  electrcsnagnetic  scattering  from  acoustic  waves,  and 
(3)  infrared  tracking  of  an  artificially  heated  bubble.  It  was  tentatively 
agreed  that  the  use  of  a  heated  bubble  as  a  sensor  did  not  merit  further 
investigation.  In  Phase  III,  therefore,  experiments  will  be  designed  only  for 
the  wind  measuring  methods  using  electromagnetic  backscattering  from  natural 
turbulence  and  from  acoustic  waves. 
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II. 


FACTUAL  DATA 


A.  Introduction 

A  definite  Arny  requirement  exists  for  remote  soundings  to  determine 
the  basic  parameters  of  the  atmosphere  without  the  need  for  probes  or  material 
sensors  located  at,  or  passing  through,  the  region  of  Interest.  The  objective 
of  this  program,  as  specified  in  Technical  Requirement  Nr  SCL-5825,  is  a  study 
of  possible  methods  for  the  determination  of  the  wind  field  of  the  atmosphere, 
directly  or  Indirectly,  by  active  or  passive  means  at  space  coordinates  remote 
from  the  location  of  the  measuring  apparatus. 

The  first  phase  of  the  work  was  a  survey  of  the  literature  and  or¬ 
ganizations  In  the  field  to  find  possible  methods  of  remote  wind  measurement. 
Nine  proposed  methods  were  found  In  the  survey,  but  a  preliminary  analysis, 
based  on  the  results  of  other  investigators,  left  only  three  of  these  methods 
for  detailed  consideration  in  Phase  II,  Analysis. 

The  analyses  in  this  phase  have  been  made  to  determine  the  feasibility 
of  using  each  of  the  three  methods  as  the  basis  for  a  remote  wind-measuring  in¬ 
strument.  Consideration  has  been  given  to  such  factors  as  physical  parameters 
Involved,  the  theory  of  operation,  available  data,  and  the  state-of-the-art  of 
radar  and  infrared  instrumentation.  Where  required  characteristics  are  not 
known,  use  has  been  made  of  mathematical  models  based  on  assumed  characteristics. 

The  three  methods  considered  are: 

1.  Scattering  from  natural  atmospheric  turbulence, 

2.  Electromagnetic  scattering  from  acoustic  waves,  and 

3.  Infrared  tracking  of  an  artificially  heated  bubble  of  air. 

The  first  method  has  previously  been  mentioned  by  veurlous  Investigators,  but  no 
detailed  examination  has  been  made.  The  second  method  has  been  demonstrated 
at  short  ranges  and  a  general  system  analysis  performed.  The  third  method  was 
suggested  by  Mr.  Harry  Moses  of  Argonne  National  Laboratory  in  a  discussion  dur¬ 
ing  the  survey  phase  and  is  not  known  to  have  been  considered  previously.  Sim¬ 
ilar  analytical  procedures  would  be  desirable  so  that  the  results  could  be  put 
on  a  more  common  basis  for  comparison.  The  variation  in  parameters  used,  avail¬ 
able  information,  and  previous  work  done,  made  it  difficult  to  follow  a  rigid 
procedure  in  the  analysis  of  each  method. 


B. 


Scattering  from  Natural  Atmospheric  Turbulence 


1.  Method:  One  of  the  moet  attractive  concepts  for  remote  measure¬ 
ment  of  wind  velocities  Involves  the  use  of  reflections  of  either  radar  or 
acoustic  signals  from  the  turbulence  which  exists  in  the  atmosphere.  A  system 
mairing  use  of  such  reflections  is  ideal  in  that  the  atmosphere  furnishes  its 
own  sensor,  in  the  form  of  turbulent  "blobs"  or  eddies,  at  the  point  where  the 
wind  velocity  is  to  be  measured.  In  principle,  the  motion  of  this  turbulent 
structure  only  need  be  observed  to  ascertain  the  motion  of  the  adjacent  atmos¬ 
phere.  The  turbulent  "blobs"  woxild  be  assumed  to  move  along  with  the  mean 
wind  speed.  Their  velocity  would  be  determined  by  DoK>ler  techniques,  because 
the  transient  nature  of  the  turbulent  structure  makes  it  unlikely  that  an  ord¬ 
inary  position-tracking  system  could  be  used  for  velocity  measurements. 

Unfortunately,  very  little  information  is  avsdlable  on  the  eu:t:ial, 
detailed  characteristics  of  low  level  turbulence.  This  is  due  partly  to  the 
lack  of  adequate  instrumentation  Etnd  partly  to  the  lack  of  a  pressing  demand 
for  such  data.  Past  Interest  in  the  subject  is  well  summeurized  by  R.  J.  Papai/ 
in  his  rather  comprehensive  survey  of  atmospheric  turbulence : 

"Anemometer  recordings  illustrate,  especially  in  the  lower 
layers  of  the  atmosphere,  tliat  the  wind  is  usually  highly  irregular. 

The  wind  velocity  is  composed  of  a  complex  of  oscillations  of  dura¬ 
tion  varying  from  a  fraction  of  a  second  to  many  minutes  and  of  an 
amplitude  >rtiich  is  often  a  substantial  fraction  of  the  average  speed. 

The  speed  of  the  wind  chauiges  not  only  from  instwt  to  instant  but 
also  from  point  to  point  of  space.  A  complete  specification  of  the 
velocity  field  over  even  a  limited  portion  of  the  atmosphere  is  im¬ 
practicable.  The  study  of  atmospheric  turb\ilence  is  primarily  con¬ 
cerned  with  mean  values,  such  as  analysis  of  the  mean  distribution 
of  momentum,  heat  and  suspended  particles  in  this  rapidly  changing 
field." 

If  it  is  assumed  that  atmospheric  turbulence  will  always  be  available 
when  needed,  there  still  remain  several  serious  obstacles  to  any  successful 
development.  In  particular,  there  are  three  basic  problems  which  must  be 
solved  before  the  concept  of  reflections  from  natural  turbulence  can  be  put 
to  practical  use  for  measuring  wind  velocities.  These  are  the  question  of 
turbulence  motion  with  the  wind,  the  effects  of  Doppler  broadening,  suid  the 
problem  of  adequate  reflection  levels  from  natural  turbulence. 

2.  The  motion  of  turbulence  with  respect  to  the  wind:  It  is  assumed 
that  the  turbulence  in  the  atmosphere  may  be  characterized  by  a  more-or-less 
random  array  of  eddies  or  "blobs."  However,  these  "blobs"  will  be  usable  for 
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■wind  ■velocity  measurements  only  if  they  are  carried  along,  or  con^vected,  with 
a  ■velocity  equal  to  the  mean  wind  ■velocity,  while  maintaining  essentially  uni¬ 
form  structure.  There  is  considerable  reason  to  doubt  that  this  should  be  so. 
The  turbulence  represents  a  random  motion  superiiqposed  on  the  mean  motion,  and 
the  turbulent  "blobs"  art  believed  to  be  continually  growing  and  dissipating, 
rather  than  remaining  in  a  fixed  condition. 

Some  experimental  work  directed  toward  a  solution  to  this  problem 
has  been  conducted  at  the  University  of  Southampton. §/  The  experiments  have 
been  directed  at  establishing  whether  the  turbulent  structure  in  a  moving  air 
stream  moves  with  the  mean  stream  velocity.  This  has  been  approached  by 
studying  the  correxatim  of  the  measured  turbulent  fluctuations  (of  velocity) 
as  a  function  of  downstream  distance  and  time.  If  the  turbulent  "blobs"  do 
beha^ve  as  quasi-permanent  bodies  \dd.cb  are  ceorrled  along  at  the  mean  stream 
velocity,  these  correlation  experiments  would  be  expected  to  show  results 
similar  to  the  curves  in  Fig.  1,  where: 

R(t)  =  correlation  between  quemtity  as  measured  at  one  station  and 
at  another  station  a  distance  x  downstream; 

T  =  time  difference  between  the  measurements  at  the  two  stations; 

X  =  distance  to  downstream  station;  and 
Vg  =  mean  stream  velocity. 

For  exan^le,  a  correlation  between  measurements  taken  at  two  stations  a  dis¬ 
tance  dj^  ap€u:t  that  peaks  at  a  ■value  of  1.0  at  time  t  =  dj/Vg  indicates 
that  the  turbulent  structure  has  remained  essentially  constant  while  being 
con^vected  the  distance  d^  with  a  velocity  equal  to  the  mean  stream  ■velocity, 
Va  . 


The  actual  experimental  results  correspond  more  closely  to  Fig.  2. 

In  these  curves,  the  correlation  coefficient  peaks  at  a  value  less  than  unity, 
8uid  the  peak  value  decreases  as  the  downstream  distance  Increases.  This  shows 
that  the  turbulent  structure  is  changing  with  time,  and  emphasizes  the  transi¬ 
tory  nature  of  the  "blob"  structure.  In  the  experiments  reported,  the  ■value 
of  the  peak  R(t)  decreased  below  0.5  in  150  psec.  or  less.  In  addition,  the 
peaks  occur  at  times  greater  than  the  corresponding  d/V  values,  indicating 
that  the  turbulent  structure  tends  to  move  more  slowly  than  the  mean  stream 
velocity. 


While  these  results  appear  to  be  discouraging  to  the  wind -measuring 
application,  several  points  concerning  the  experimental  work  should  be  men¬ 
tioned.  First,  the  experiments  were  performed  downstream  from  a  Jet  of  air 
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Fig.  1  -  Results  Expected  from  Idealized  Correlation  Experiment 


Fig.  2  -  Results  of  Actual  Correlation  Experiment  in  a  Moving  Air  Stream 

(after  Richards  and  Williams!/) 


expelled  through  a  1-ln.  orifice,  cmd  thus  it  may  be  questioned  whether  the 
results  are  applicable  to  turbulence  in  the  free  atmosphere.  Secondly,  the 
correlation  measiurenents  were  performed  on  velocity  fluctuations,  whereas  it 
is  the  behavior  of  the  index  of  refraction  fluctuations  that  is  of  Interest 
in  the  electromagnetic  reflection  problem  and  that  of  the  acoustic  Impedance 
in  the  acoustic  case. 

Consequently,  it  appears  necessaury  to  perform  similar  correlation 
experiments  in  the  free  atmosphere  using  measurements  of  refractive  index  or 
acoustic  impedance  before  an  answer  to  the  problem  of  convection  of  turbulence 
can  be  given.  If  a  similar  decrease  of  the  peak  correlation  coefficient  is 
observed,  there  may  still  be  a  possibility  of  using  this  system,  providing  the 
peak  R(t)  remains  close  to  unity  for  a  time  t  long  enough  to  permit  the 
required  Doppler  frequency  measurement.  If  this  happens,  the  subsequent  decay 
of  R(t)  would  not  affect  the  operation  of  the  proposed  wlnd>measurlng  system. 

3.  The  effect  of  Doppler  broadening:  Another  obstacle  to  the  use 
of  Doppler  measureinents  of  turbulence  motion  with  the  wind  is  the  effect  of 
the  random  motion  associated  with  the  turbulence.  In  sux  ideall:’.ed  situation, 
the  turbulent  "blobs"  could  be  regarded  as  well-defined  bodies  moving  at  a 
viniform  velocity  with  the  wind.  The  spectrum  of  Doppler  frequencies  in  a 
signal  reflected  from  the  moving  "blobs"  would  then  consist  of  a  single  line 
at  the  frequency  corresponding  to  the  wind  component  along  the  line  of  sight. 

In  the  actual  situation,  there  is  a  random  turbulent  motion  superiaposed  on 
the  mean  wind  speed.  The  resulting  spectrum  of  velocities  will  give  rise  to 
a  wide  spectrum  of  Doppler  frequencies.  This  effect  is  illustrated  in  Fig.  3; 
a  similar  phenomenon  is  the  well  known  Doppler  broadening  encountered  in 
spectroscopy,  >diere  the  random  thermal  motions  of  the  molecules  give  rise  to 
a  broadening  of  the  spectral  lines. 

The  Inoportance  of  this  Doppler  broadening  effect  is  difficult  to 
evaluate  without  some  experimental  evidence.  A  set  of  measurements  performed 
by  Emerson  Research  Laboratories^/  in  which  the  Doppler  frequencies  in  an 
acoustic  signal  reflected  from  atmospheric  turbulence  were  determined  provides 
some  evidence  of  this  effect.  For  wind  speeds  up  to  25  knots,  a  Doppler  spec¬ 
trum  width  corresponding  to  +10  knots  was  not  uncommon.  While  it  is  not  cer¬ 
tain  that  all  of  this  error  can  be  attributed  to  Doppler  broadening,  these 
experiments  do  at  least  provide  a  clue  to  the  possible  importance  of  this 
effect. 


4.  Ihe  problem  of  adequate  reflection  levels  from  turbulence:  The 
third  problem  is  that  of  obtaining  from  atmospheric  turbulence  a  reflected  sig¬ 
nal  which  is  strong  enough  to  be  detected.  Although  the  work  by  Emerson  demon¬ 
strated  the  possibility  of  detecting  acoustic  reflections  at  short  ranges  (72.5 
feet)  the  experimental  configuration  was  not  "remote,"  but  resembled  a  "forward- 
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Fig.  3  -  Effect  of  Turbulent  Motion  on  the  Doppler  Frequency  Spectrum 


scatter"  propagation  path.  It  Is  believed  that  the  acoustic  method  vlU  be 
limited  to  very  short  ranges,  because  of  background  noise  at  low  frequencies 
and  atmospheric  absorption  at  high  frequencies. 

There  Is  considerable  evidence  that  radar  reflections  from  turbulence 
have  been  detected  by  existing  radar  systems.  There  are,  for  example,  frequent 
reports  of  radar  "angels"  which  can  be  attributed  only  to  turbulent  Inhomogene- 
Itles  of  the  atmosphere.  There  are  reported  Insteuices  In  which  such  "angels" 
have  been  tracked  by  radar  for  periods  of  several  minutes.  Further  evidence 
In  support  of  this  concept  may  be  found  In  the  theories  of  tropospheric  scatter 
propagation,  which  are  usually  based  on  a  mechanism  of  scattering  from  an  ele> 
vated  turbulent  layer.  Radar  reflections  from  periodic  dlst\irb8mces  In  the 
atmosphere  (see  Section  II-C)  sure  greatly  erihanced  at  certsdn  ratios  of  radar 
to  acoustic  wavelengths,  leading  to  the  belief  that  radar  wavelengths  might 
be  varied  to  provide  maximum  returns  for  a  given  scale  of  turbulence. 

None  of  this  evidence  may  be  regarded  as  establishing  the  ability  of 
radar  to  detect  atmospheric  turbulence  consistently  (in  the  absence  of  precipi¬ 
tation  returns).  However,  the  possibility  is  so  attractive  that  both  the 
Weather  Bureau^  and  the  Vteather  Radar  Group  under  Dr.  David  Atlas  at  the  Air 
Force  Cambridge  Research  Laboratories^/  have  recently  initiated  projects  to 
investigate  this  possibility  further.  This  is  also  an  area  of  great  concern 
in  connection  with  aircraft  warning  for  avoidance  of  clear-air  turbulence. 
Considerable  progress  may  be  expected  toward  a  satisfactory  answer  to  the 
question  of  detectability  of  atmospheric  turbulence  in  the  near  future. 

One  of  the  most  important  limitations  on  the  use  of  atmospheric 
turbulence  for  wind  velocity  measurements  will  be  the  attainable  range.  At 
present  this  is  difficult  to  estimate,  since  no  satisfactory  theory  of  radar 
reflections  from  atmospheric  turbulence  exists.  Perhaps  the  best  estimate 
can  be  obtained  from  the  observations!/  of  "angels",  which  in  most  cases  were 
within  5,000  yards  of  the  radar.  The  maximum  range  to  >rtiich  atmospheric  txir- 
bulence  could  be  consistently  detected  under  a  veuriety  of  meteorological  con¬ 
ditions  must  be  regaurded  as  at  present  unknown. 

5.  Conclusions:  The  concept  of  wind  measurement  by  radar  determina¬ 
tion  of  mean  turbulence  velocity  is  quite  attractive  due  to  the  use  of  a 
natural  phenomenon  as  the  remote  sensor.  Considerably  more  detailed  informa¬ 
tion  will  be  required,  however,  regarding  the  8w:t\ial  characteristics  of  low- 
level  turbulence,  if  the  feasibility  of  the  concept  is  to  be  established 
analytically.  In  particular,  data  will  be  needed  on: 

(1)  Whether  turbulence  actixally  moves  with  the  wind; 

(2)  The  variations  of  scale  for  turbulence  in  the  first  mile  of 

the  atmosphere; 
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(3)  The  range  of  index  of  refraction  variations  to  be  expected; 

(4)  The  range  of  random,  internal  velocities  of  the  turbulent 

cells;  and 

(5)  The  probability  of  having  turbulence  available  under  a  vide 

range  of  atmospheric  conditions. 

Successful  application  of  this  concept  would  then  require  that  the 
radar  state-of-the-art  be  advanced  to  the  point  where  the  reliability  of  back- 
scatter  measurements  approached  that  of  present  forward- scatter  conanunication 
systems.  If  this  point  is  reached,  it  is  believed  that  statistical  processing 
of  the  data  will  enable  the  resolution  of  useful,  if  not  exact,  velocity  in¬ 
formation  from  the  ambiguities  due  to  Doppler  broadening  and  the  transient 
buildup  and  decay  of  the  turbulent  cells. 


C.  Electromagnetic  Scattering  from  Acoustic  Waves 

1.  Method:  A  technique  using  electromagnetic  scattering  from 
acoustic  waves  for  the  remote  measurement  of  low-level  wind  velocity  and  tem¬ 
perature  profiles  has  previously  been  developed  and  reported!/  by  Midwest 
Research  Institute,  Through  this  technique,  air  velocity  is  determined  by 
measuring  the  Doppler  frequency  shift  in  an  electromagnetic  wave  that  has  been 
reflected  from  an  acoustic  wave  propagating  in  the  atmosphere.  Ibis  Doppler 
frequency  corresponds  to  a  velocity  that  is  the  sum  of  the  local  sound  velocity 
and  the  coiqponent  of  wind  velocity  in  the  direction  of  acoustic  propeigation. 
Since  the  velocity  of  sound  is  a  function  of  air  temperature,  the  Doppler  fre¬ 
quency  measurement  contains  both  wind  and  temperature  information. 

The  electromagnetic-acoustic  (EMAC)  probe  consists  of  a  high-intensity 
directional  acoustic  source  and  a  Doppler  radar  system  as  shown  in  Fig.  4.  To 
obtain  complete  wind  and  temperature  information,  the  volume  of  interest  must 
be  probed  from  different  directions;  ordinarily,  this  will  require  a  separate 
EMAC  probe  for  each  sounding  direction.  In  the  most  general  case,  four  sound¬ 
ing  directions  (hence  four  probes)  will  be  required  to  determine  four  quanti¬ 
ties;  three  vector  components  of  the  wind  velocity  and  temperature.  If  cer¬ 
tain  assumptions  can  be  made,  the  number  of  probes  can  be  reduced.  These  are: 
if  the  vertical  component  of  the  wind  velocity  is  negligible,  then  the  number 
of  probes  needed  is  reduced  to  three;  and  if  the  air  temperature  is  constant 
throughout  the  volume  of  interest,  another  probe  may  be  eliminated  in  favor  of 
a  more  economical  meems  of  temperature  measurement. 
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In  cases  where  the  vertical  wind  velocity  Is  small,  a  single  EMAC 
prohe  can  he  used  to  measure  temperature  profiles  only,  along  a  vertical  path. 
Thus,  a  single  probe  could  be  used  for  sounding  low-level  temperature  Inver¬ 
sions  in  the  atmosphere.  An  accuracy  of  1  to  2”F  and  a  height  of  1,500  ft. 
appeea*  to  be  reasonable  performance  criteria  for  this  application. 

In  the  reported  experimental  system,  the  acoustic  signal  treuismitted 
is  a  e-ms.  pulse  of  22-kc.  frequency.  The  acoustic  intensity  is  about  140  db. 
(re  0.0002  dyne/cm^)  at  10  ft.  from  the  source,  and  the  beamwldth  is  about 
5  degrees.  The  radar  system  is  a  3-w.  CW  Doppler  radar  operating  at  10,000  me.; 
receiver  sensitivity  is  -96  dbm.  The  acoustic  wavelength  is  1.5  cm.,  while 
that  of  the  radar  signal  is  3  cm.;  this  one-to-two  wavelength  relationship 
must  be  maintained  to  obtain  usable  reflections. 

With  this  experimental  system,  radar  echoes  from  the  acoustic  pulse 
have  been  observed  out  to  a  distance  of  about  90  ft.  The  range  is  limited 
chiefly  by  atmospheric  absorption  of  the  acoustic  energy  \rtiich,  at  22  kc., 
exceeds  0.2  db/ft.  By  operating  the  system  at  lower  frequencies  where  absorp¬ 
tion  is  less,  the  range  may  be  greatly  extended.  For  Instance,  at  5.5  kc. 

(6  cm.),  the  absorption  is  only  0.06  db/ft  and  ranges  up  to  1,500  ft.  should 
be  attainable. 

A  detailed  analysis  of  the  EMAC  probe  wind  measuring  system  has  been 
made§/  and  is  included  in  this  report  as  Appendix  A.  Consequently,  enough  in¬ 
formation  about  the  probe  is  available  to  permit  good  estimates  of  such  per¬ 
formance  criteria  as  range  of  measurement,  accuracy  of  measurement  in  speed 
and  direction,  resolution,  response  time,  complexity  and  reliability.  Dis¬ 
cussions  on  these  performance  criteria  are  given  in  the  following  sections. 

2.  Maximum  range  of  measurement:  The  msucimum  operating  range  of  a 
single  EMAC  probe  is  determined  by  the  frequencies  employed,  transmitter  power, 
receiver  sensitivity,  condition  of  the  atmosphere,  acoustic  pulse  length  Md 
intensity,  and  the  alignment  of  radar  and  acoustic  wavefronts.  The  analysis  has 
led  to  the  conclusion  that  for  practical  choices  of  pareuneters  a  maximum  use¬ 
ful  range  of  1,500  ft.  would  be  expected  for  a  single  probe  operating  at  2,500 
me.  Recent  e:q)erimental  work®/  however  indicates  that  operation  may  be  pos¬ 
sible  at  frequencies  as  low  as  300  me.,  with  attendant  Increases  in  range. 

The  multi-probe  arrays  required  for  determining  complete  wind-vector  data  will 
be  further  limited  in  range  because  of  scanning  requirements;  the  chief  limit¬ 
ing  factor  here  is  the  loss  of  accuracy  of  the  multi-probe  arrays  at  the  longer 
ranges.  This  is  due  to  geometric  dilution  at  the  small  intersection  angles. 
Examination  of  the  results  of  the  error  analysis  (Appendix  A)  shows  that  the 
maximum  accurate  range  of  a  multi -probe  array  will  be  comparable  to  the  base¬ 
line  length  available  for  siting  the  several  probes,  amd  at  most  will  be 
somewhat  less  than  the  maximum  range  of  the  individual  probes. 
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5*  Accuracy  of  vind  measurement;  The  probe  arrays  suitable  for 
determining  winds  remotely  have  been  designed  to  provide  as  outputs  the  values 
of  two  horizontal  wind  components,  which  may  be  regaz^d  as  the  range  and 
cross-wind  components,  respectively.  The  typical  accuracy  of  the  measured 
wind  components  will  be  about  ±  0.5  nph  to  t  1.0  mph;  the  exact  values  will 
vary  along  the  selected  profile  as  the  antenna  sighting  angles  vary.  Combin¬ 
ing  these  components  into  the  resultant  wind  vector  results  in  acc\u‘acles  for 
the  total  wind  speed  of  t  1.0  to  t  1.5  iqph. 

The  accuracy  of  determination  of  the  wind  direction  is  obviously 
limited  by  the  accuracy  with  \diich  the  wind  vector  can  be  ireconstructed  from 
the  measured  components.  If  the  actual  wind  speed  is  only  1  nqph,  then  the 
errors  in  the  measured  components  will  mske  accurate  determination  of  the  wind 
direction  almost  Impossible.  However,  as  the  wind  speed  increases  the  accuracy 
of  direction  determination  eJ.80  increases.  For  example,  referring  to 
Fig.  5,  the  azimuth  error  at  a  wind  speed  of  10  mph  would  be  about  t  6*.  This 
does  not  include  siting  errors  and  azimuth  errors  in  the  alignment  of  the 
antenna  itself. 

Three-dimensiCTial  resolution  of  measurement  position;  For  the 
BMAC  probe,  the  spatial  resolution  of  the  wind  velocity  measurements  is  a 
function  of  many  system  design  parameters.  The  resolution  of  a  single  probe, 
Fig.  6a,  in  the  radial  direction  is  determined  by  the  spatial  length  of  the 
acoustic  pulse  and  by  the  distance  the  acoustic  pulse  travelled  during  the 
Doppler  frequency  measurement.  An  increase  in  the  length  of  the  acoustic 
pulse  is  desirable  because  it  Increases  the  radar  reflection  coefficient  of 
the  acoustic  wave  packet;  on  the  other  hand  a  shorter  pulse  provides  greater 
spatial  resolution.  The  analysis  Indicated  that  a  satisfactory  compromise  can 
be  obtained  by  using  an  acoustic  pulse  length  of  about  20  to  25  ms.  which 
results  in  a  range  resolution  of  about  25  ft. 

In  the  two  dimensions  perpendicular  to  the  radial  direction,  the 
resolution  is  governed  primarily  by  the  width  of  the  acoustic  and  reuiar  beams. 
For  a  fixed  angular  beamwldth  the  resolution  distance  will  increase  in  propor¬ 
tion  to  the  range.  For  exanple,  at  a  remge  of  250  feet  a  conical  acoustic 
beamwidth  of  6*  provides  a  resolution  of  about  25  ft.,  while  at  a  range  of 
1,000  ft.  the  resolution  would  be  about  100  ft.  In  principle,  the  beamwldth 
can  be  reduced  to  Increase  the  resolution,  but  problems  of  increasing  antenna 
size  and  the  difficulty  of  keeping  the  acoustic  and  radar  beams  aligned  will 
be  limiting  factors. 

When  a  multi-probe  array  is  used  to  obtain  complete  wind  vector  in¬ 
formation,  some  loss  of  resolution  will  result  because  the  probes  must  examine 
the  volume  of  interest  from  different  directions. 
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IF  AVx  *  0.7  m 
AVy  =  0.7  MM 

THEN  AV  =  1.0  MM 
FOR  V--  10. 0  MM 

^  V  10.0  MM  O-'WDIMI 

Fis»  5  -  Typical  Azimuth  Error  in  Wind  Direction  for  Wind  Speed  of  10  MPH 
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Figure  6b  illustrates  a  typical  situation  for  a  three-probe  array;  it  is  evi¬ 
dent  that  the  resolution  obtainable  with  a  multi-probe  array  will  be  comparable 
to  the  largest  dimension  of  any  of  the  individual  scanning  beams  at  the  loca¬ 
tion  being  examined. 

5.  Measurement  and  display  response  time;  The  EM/VC  probe  should 
actually  be  cheuracterlzed  by  tvo  response  times.  The  first  is  the  time  re¬ 
quired  to  make  a  measurement  of  the  Doppler  frequency  corresi>onding  to  the 
velocity  of  the  acoustic  vave  packet  at  any  given  location,  assuming  that  the 
acoustic  vave  packet  has  already  reached  this  location.  This  time  vlll  be  of 
the  order  of  a  few  milliseconds,  depending  on  the  pairtlcular  method  of  fre¬ 
quency  measurement  employed  and  the  required  siccuracy.  This  response  time  af¬ 
fords  a  measure  of  the  ability  of  the  EMAC  probe  to  respond  to  rapid  fluctua¬ 
tions  in  air  velocity  or  fluctuations  of  temperature,  if  provision  is  made  for 
extreuitlng  temperature  information  from  the  measurement  of  sound  velocity.  The 
second  response  time  of  Interest  is  that  required  for  the  scanning  of  a  given 
profile  or  a  volume  of  space,  uhich  will  be  determined  primarily  by  the  size 

of  the  volume  of  interest  and  by  the  propagation  time  of  the  acoustic  vave. 

Ihe  acoustic  signal  has  a  velocity  of  1,100  ft/sec,  so  the  time  required  to 
scan  a  radial  path  vith  one  EMAC  probe,  out  to  a  distance  of  1,100  ft.,  vould 
be  1  sec.  Examination  of  a  larger  volume,  \diere  several  scans  in  different 
directions  are  required,  or  scanning  a  single  profile  vith  several  probes,  vlll 
require  a  correspondingly  longer  time. 

As  an  illustration  of  this  last  point,  a  detailed  scanning  program 
suitable  for  one  particular  application  has  previously  been  vorked  out.iS/ 

The  values  of  the  vlnd  vector  at  five  points  out  to  a  maximum  range  of  300  yd. 
were  determined  by  a  three-probe  array.  The  geometry  of  the  assumed  situation 
is  shovn  in  Fig.  7.  The  required  scanning  time  vas  found  to  be  about  5  sec. 
for  this  particular  case. 

6.  Data  presentation  form;  The  experimental  versions  of  the  EMAC 
probe  which  have  been  used  to  date  have  had  only  simple  forms  of  data  display. 
These  have  consisted  of  oscilloscope  presentation  of  the  Doppler  signal,  and 

a  gated  counting  circuit  used  for  measurement  of  the  Doppler  frequency.  In 
operational  versions  of  the  EMAC  probe,  the  display  vould  present  either  range 
and  crossvlnd  coinponents,  or  resultant  wind  speed  and  direction,  or  both.  The 
display  vould  probably  be  in  digital  form,e8pecially  if  the  equations  \rtiich 
must  be  solved  to  obtain  the  vind  velocity  components  from  the  measured  Doppler 
frequencies  are  solved  by  a  small  digital  computer. 

Other  data  which  could  be  displayed  Include  air  temperature  (actually 
the  so-called  "virtual  temperature"  obtained  from  the  measured  sound  velocity) 
and  a  weighted  wind-correction  factor  suitable  for  use  in  fire  control-  or 
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Fie.  7  -  A  Typical  Wind  Profile  Problem 
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rocket  launching  applications.  Where  the  wind-correction  factor  is  desired, 
it  could  be  calculated  by  the  same  cotqputer  used  to  calculate  the  wind 
cooponents. 

In  cases  where  more  cooplete  wind  field  data  are  desired,  problems 
of  data  storage  and  presentation  will  be  encountered.  These  problems  axe  not 
unique  to  the  EMA.C  probe,  but  would  be  encountered  in  any  system  intended  for 
measuring  a  detailed  wind  field.  The  data  logging  equipment  at  the  Cedar  Hill 
tower, ^whlch  records  wind  and  tenperature  Information  at  only  12  points, 
gives  some  idea  of  the  scope  of  the  problem. 

7.  Complexity  of  the  required  instrumentation:  A  block  diagram  of 
a  proposed  EMAC  probe  system  is  shown  in  Figs.  8  and  9.  Although  the  conplex- 
ity  of  the  instrumentation  is  difficult  to  estimate  from  block  diagrams,  cer¬ 
tain  conclusions  can  be  drawn.  The  complexity  of  the  radar  portion  of  the 
system  will  be  comparable  to  that  of  the  battlefield  surveillance  radars  now 
in  existence.  The  acoustic  system  presents  added  complexity;  however,  in  the 
final  system  design,  a  fixed  acoustic  frequency  will  be  used,  requiring  a  very 
simple  acoustic  source. 

8.  Capability  for  total  wind  field  display;  A  single  EMAC  probe  can 
be  rotated  in  azimuth  to  examine  any  location  within  the  range  of  the  probe, 

up  to  elevation  angles  of  approximately  IS  degrees  (where  the  error  introduced 
by  neglecting  the  vertical  motion  may  become  significant).  However,  the  single 
scanning  probe  gives  accurate  measurements  only  when  the  wind  is  stratified  so 
that  the  wind  velocity  may  be  considered  to  be  a  function  of  height  only. 

Within  the  limits  imposed  by  reuige  and  elevation  coverage,  the  single-probe 
system  can  provide*  deteii led  Information  about  the  variation  of  the  v  nd  veloc¬ 
ity  with  height. 

The  ability  of  the  EMAC  probe  system  to  display  a  "simultaneous" 
wind  field  is  governed  by  the  response  time  considerations  previously  dis¬ 
cussed.  Assuming  that  the  volume  of  interest  is  to  be  scanned  in  a  sequential 
fashion,  a  time  lapse  between  measurements  at  various  points  in  the  volume 
will  be  inevitable. 

9.  Reliability;  Since  only  experimental  models  of  the  EM\C  probe 
have  been  built,  there  is  little  experience  upon  \dilch  to  base  Judgment  of  its 
expected  reliability.  The  radar  portion  of  the  probe  should  coopare  in  relia¬ 
bility  to  small  Doppler  radars  now  in  existence.  The  reliability  of  the  siren 
acoustic  source  may  present  difficulties;  moving  parts  consist  of  a  high-speed 
siren  and  a  magnetically-operated  valve  used  to  pulse  the  compressed  8d.r  to 
the  siren.  The  development  of  an  acoustic  source  requiring  no  moving  peurts 
would  be  of  considerable  advantage. 


17 


18  - 


Block  Diagram  of  Proposed  IMAC  Pulse  Doppler  Radar  System 
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Fig.  9  -  Block  Diagram  of  EMAC  Probe  Acoustic  Source 


10.  Conclusions;  The  EMAC  probe  at  microwave  frequencies  appears 
to  be  a  useful  technique  for  remote  wind  measurements  at  remges  up  to  approxl- 
mately  1,500  ft.  \dien  operating  at  2,500  me.  Its  extremely  fast  response  to 
velocity  variations  and  its  sensitivity  at  low  velocities  make  it  quite  useful 
for  measurement  of  mlcrometeorologlcal  parameters  such  as  gusts,  turbulence, 
and  the  detailed  thermal  structure  of  the  atmosphere. 

The  resolution  of  the  EMAC  probe  In  time  and  space  depends  In  a  eexn- 
pUcated  way  on  several  variables  such  as  beamwldths,  sensitivity  of  available 
receivers,  range,  and  number  of  probes  needed.  A  single  probe  has  a  spatial 
resolution  of  25  ft.  or  less  and  a  response  time  of  about  50  ms.  Vlhen  an  array 
of  probes  Is  used  to  obtain  coiqplete  wind  axid  temperature  data,  the  time  re¬ 
quired  to  collect  and  Interpret  the  data  from  the  various  soundings  will  de¬ 
grade  the  response  time  by  an  order  of  loagnitude.  In  addition,  some  loss  of 
spatial  resolution  will  result  if  the  volume  of  intersection  of  the  several 
soundings  has  dimensiois  greater  than  25  ft.  The  expected  accureicies  In  the 
measurements  will  also  depend  on  the  particular  design,  but  values  of  t  1  mph 
In  wind  velocity  and  t  1*F  in  tenperature  are  typical. 

Its  most  severe  limitation  Is  the  meocimum  range  quoted  above.  This 
figure  Is  an  extension  of  experimental  data  at  higher  frequencies  and  is  be¬ 
lieved  to  be  the  maxliiium  ssUTe  extrapolation  of  probe  characteristics  without 
additional  experimental  verification.  The  analysis,  yhttn  carried  to  even 
lower  system  frequencies,  predicts  greater  usable  ranges  due  to  reduced 
acoustic  attenuation,  and  recent  e;q>eriments  indicate  that  operation  may  be 
possible  at  frequencies  as  low  as  300  me.  However,  additional  unknown  factors 
must  be  investigated  before  meaningful  analytical  extrapolations  can  be  made. 
These  factors  include  attenuation  due  to  turbulent  scattering  of  the  acoustic 
signed,  random  scattering  from  turbulence  "noise"  of  the  same  scale  as  the 
acoustic  wavelength,  and  the  effects  of  acoustic  beam  displacement  by  the  wind 
at  Icxig  ranges.  Here,  as  in  the  preceding  method,  an  increased  knowledge  of 
low-level  turbulence  structure  is  needed  before  accurate  calculations  can  be 
made  of  expected  system  performeuice . 

Although  this  method  is  the  only  one  of  the  three  that  has  been 
demonstrated  experimentally,  considerable  development  is  needed  before  a  com¬ 
plete  long-range  wind-field  measuring  instrument  is  produced.  Further  ejqperi- 
mental  work  should  be  performed  to  verify  the  predicted  maximum  usable  ranges 
for  lower  frequencies  euid  determine  the  optimum  system  parameters  for  maximum 
range.  In  addition,  an  experimental  system  could  now  provide  a  more  precise 
verification  of  the  EMAC  probe  theory,  permit  precise  measurements  of  reflec¬ 
tion  coefficient  and  scattering  cross-section  of  the  acoustic  waves,  and  pro¬ 
vide  for  the  first  time  an  instrument  capable  of  making  area,  rather  than 
point, measurements  of  low-level  turbulence  and  wind  profiles. 
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D.  Infrared  Tracking  of  an  Artificially  Hteated  Bubble 


1.  Method;  This  method  would  use  remote  means,  such  as  Intersecting 
acoustic,  Infraz^d,  or  microwave  beams  which  could  be  focused  on,  and  absorbed 
by  a  remote  volume  of  air,  to  raise  Its  tenqperature  above  ambient.  The  heated 
volume  of  air,  referred  to  as  a  bubble,  would  rise  in  the  atmos^ere  at  a  rate 
dependent  on  Its  temperature  difference  and  the  local  lapse  rate.  As  it  rises. 
Its  path  would  be  Influenced  by  the  locad  wind  field  In  the  same  manner  as  a 
balloon.  Theoretically,  the  bubble  could  be  tracked  by  a  multiple  Infrared 
theodolite  system  to  obtedn  its  position  as  a  function  of  time,  from  which  the 

\dnd  profile  could  be  determined. 

2.  Generation  of  bubble;  The  proposed  method  of  bubble  generation 
Involved  Intersecting  beams  of  some  form  of  propagating  energy  -  acoustic, 
optical,  or  microwave  -  capable  of  being  focused  on  the  region  to  be  heated. 

Such  energy  transfer  and  concentration  techniques  have  been  used  with  success 
vhere  the  distances  from  the  transducers  to  the  Intersection  region  are  com¬ 
parable  to  the  spacing  between  transducers,  and  to  the  dimensions  of  the  trans¬ 
ducers  as  shown  In  Fig.  10.  Such  a  configuration  is  obviously  lnq>ractlcal  If 
ranges  of  miles  or  even  hundreds  of  yards  are  being  considered,  due  to  the  re¬ 
quired  size  and  spacing  of  radiating  elements  or  transducers.  A  more  xisable 
configuration  as  shown  in  Pig.  11  would  permit  realistic  spacing  of  trans¬ 
ducers,  but  would  present  additional  problems  in  accurate  focusing  and  aiming 
to  achieve  the  desired  Intersection,  in  poor  definition  of  the  intersection 
point  In  the  radial  direction,  and  in  obtaining  the  reqxdred  rates  of  energy 
input  to  the  bubble. 

For  the  geometry  of  Pig.  11,  if  it  is  assumed,  ideally  that  the  beams 
axe  perfectly  collimated,  that  they  intersect  in  a  common  volume,  and  that  no 
power-sensitive,  nonlinear  absorption  phenomenon,  such  as  ionization,  are  in¬ 
volved,  the  power  input  to  the  bubble  is  given  by 

.  od  ,  .ill 

Pin  -  nPtlO  10  ( 1-10  ^  )  ,  (1) 


>diere  n  =  the  number  of  transducers, 

P.^  =  the  power  transmitted  by  each  transducer, 

0  =  the  attenuation  of  the  radiation  in  the  atmosphere  in  db/unit 
length. 
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Fig.  11  -  Long  Range  Intersection  Geometry 
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d  >  the  distance  to  the  intersection^  and 

t  >  the  length  of  the  intersection  volume. 

Since  the  absorption  is  linear  in  decibels ^  emd  d  is  much  greater 
than  t  ,  much  more  of  the  available  power  will  be  absorbed  in  the  path  between 
bubble  and  transducer j  than  in  the  bubble.  As  an  example,  if  microwave  energy 
at  24  Gc  is  used,  a  ■  0.2  db/km  for  adr  with  1  per  cent  water  vapor.  For  a 
10  meter  bubble  at  1  km.,  4.5  per  cent  of  the  power  would  be  absorbed  in  the 
path  emd  0.043  per  cent  in  the  bubble.  If  a  frequency  such  as  190  Gc,  with  a 
theoretical  a  >  30  db/km,  is  chosen  to  get  more  absorption  in  the  bubble,  the 
power  loss  in  the  path  will  be  99.9  per  cent  with  only  0.00667  per  cent  ab¬ 
sorbed  in  the  bubble.  The  intervening  paths  would  obviously  be  heated  much 
more  than  the  bubble,  even  if  n  becooes  quite  large. 

Considering  the  highly  optimistic  assumptions  concerning  this  method, 
it  is  believed  that  in  the  actual  case  it  would  be  completely  impractical. 
Despite  the  impracticability  of  this  method  of  generation,  euialyses  were  made 
of  the  bubble  behavior  and  detection  processes,  since  it  is  believed  that  the 
information  will  be  of  value  in  situations  where  other  bubble  generation 
schemes  may  be  usable. 

3.  Analysis  of  a  heated  bubble  of  gas  in  the  atmosphere;  The  pur¬ 
pose  of  this  analysis  is  to  determine  the  feasibility  of  obtaining  wind  velocity 
measurements  by  passive  infrared  tracking  of  a  gas  bubble  which  is  heated  6uid 
released  in  the  atmosphere.  The  feasibility  depends  primarily  on  the  ability 
of  the  bubble  to  maintain  sufficient  radiance, as  it  rises  through  the  atmos¬ 
phere,  to  be  detectable  at  ground  level. 

The  analyses  of  the  thermodynamic  behavior  of  the  bubble  and  the 
infrared  detection  system  are  detailed  in  Appendix  B  and  discussed  briefly  in 
the  following  paragraphs. 

The  model  used  is  a  simple  one.  It  is  assumed  that  the  bubble  rises 
through  the  atmosphere  adiabatically  and  that  the  atmosphere  does  not  mix  with 
the  bubble  or  exert  drag  or  shear  forces  on  it.  A  consequence  of  these  assump¬ 
tions  is  that  the  model  bubble  will  radiate  more  power  than  would  be  radiated 
by  a  bubble  in  the  true  atmosphere. 

The  analysis  was  made  for  assumed  atmospheres  of  three  different 
lapse  rates;  isothermal,  l*C/l00  meters,  and  2*C/l00  meters.  For  reasons 
pointed  out  in  Appendix  B,  a  favorable  choice  of  wavelength  for  detection  is 
4.1  microns. 
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A  30  gm.  bubble  at  a  height  of  1;500  meters,  having  a  tengperature  of 
200*K  above  the  surrounding  air,  will  produce  a  radiance  of  about  10'^  watts/ 
meter^ -micron  at  the  earth's  surface  at  4.1  microns.  Ihe  background  noise  at 
the  earth's  surface  is  estimated  as  between  0.1  and  2.0  watts/meter^-steradlan- 
micron.  For  a  ground-based  detector  system  having  a  beamvidth  of  12  min.  of 
6orc,  the  solid  angle  of  view  is  9.5  x  10*^  stereidlans.  Thus,  the  maximum 
estimated  noise  is  1.9  x  10'^  watts/meter^-mlcron,  that  is,  the  same  order  of 
magnitude  as  the  signal.  On  the  basis  of  the  model  used,  the  bubble  tempera¬ 
ture  at  ground  level  would  only  be  a  few  degrees  Ifelvln  above  its  temperature 
at  1,500  meters. 

4.  Conclusions ;  An  analysis  of  bubble  generation  indicates  that  the 
original  proposal  for  heating  a  volume  of  air  with  intersecting  beams  of  energy 
does  not  appear  feasible  at  distances  ^ich  are  large  in  terms  of  transducer 
size  and  spacing.  This  eliminates  the  "remote"  characteristics  of  the  system, 
as  applied  to  this  study.  Although  the  bubble  could  be  generated  by  other 
means,  such  as  a  distant  shell  burst  or  short  duration  flare,  it  is  believed 
that  presently  available  techniq.'ues  such  as  tracked  pilot  balloons  would,  offer 
sinqpler,  more  reliable  operation. 

A  parallel  analysis  of  bubble  characteristics  was  made, assuming  that 
the  bubble  was  available  in  the  atmosphere.  The  thermodynamic  model  used  to 
represent  a  hot  bubble  in  the  analysis  was  greatly  simplified  to  make  possible 
a  reasonable  mathematical  treatment.  The  simplifications,  in  general,  led 
toward  an  optimistic  view  of  the  bubble's  lifetime  and  behavior.  The  simple 
model  shows  that  the  atssumed  bubble  could  retain  its  identity  as  a  detectable 
sensor  to  a  height  of  1,500  meters.  Additional  paurameters  such  as  radiation 
loss,  diffusion,  turbulent  mixing,  and  shear  and  drag  forces  could  be  included 
in  a  much  more  complex  analysis,  but  the  resulting  accuracy  would  suffer  from 
the  interaction  and  uncertainties  of  the  additional  model  assunptlons.  It  is 
believed  that  these  parameters  will,  in  general,  degrade  the  bubble  configura¬ 
tion,  causing  it  to  lose  its  identity  in  the  first  few  meters  of  its  travel. 
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Ill . 


CONCLUSIONS 


The  following  conclusions  are  based  on  the  results  of  the  analyses 
made  during  this  phase  on  the  three  remote  wind  measurement  methods. 


A.  Scattering  From  Natural  Turbulence 

1.  Sufficient  Information  on  the  actual  characteristics  of  low- 
level  turbulence  was  not  available  for  an  adequate  analysis  to  be  made  of  the 
feasibility  of  the  concept  of  back-scattering  from  turbulence . 

2.  Investigations  of  this  method  should  be  continued  because  of 
the  advantages  of  using  a  natural  phenomenon  as  a  se'isor. 

3.  For  purposes  of  further  analysis,  experimental  data  should  be 
obtained  regarding: 


a.  The  motion  of  turbulence  relative  to  the  wind, 

b.  The  dimensional  scale  of  low-level  turbulence, 

c.  The  refractive  Index  variations  in  low-level  turbulence, 

d.  The  random  turbulent  velocities  and  their  contribution 
to  Doppler  broadening,  and 

e.  The  probability  of  occurrence  of  detectable  turbulence. 

4.  Present  radar  emd  acoustic  systems  lack  the  required  sensitivity 
for  use  with  this  concept. 


B.  Electromagnetic  Scattering  From  Aroustlc  Waves 

1.  This  method  is  feasible  for  the  remote  measurement  of  wind  at 
short  ranges. 


2.  Experimental  data  are  needed  to  verify  existing  predictions  of 
capabilities  before  further  theoretical  extrapolations  can  be  reliably  accom¬ 
plished  . 


3 .  Further  improvements  in  instrumentation  are  needed  to  achieve 
the  predicted  capabilities  and  provide  data  to  support  existing  theory. 
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4.  Experimental  investigations  should  be  made  to  determine  how 
maximum  range  is  affected  by  such  factors  as: 


a.  Attenuation  due  to  turbulent  scattering  of  the  acoustic 

wave, 

b.  Radar  scattering  frcm  turbulence  of  the  same  scale 
as  the  acoustic  wavelength,  and 

c.  Acoustic  beam  displacement  by  the  wind. 


C.  Infrared  Tracking  of  an  Artificially  Heated  Bubble 

The  generation  of  a  well  defined  heated  region  by  intersecting 
beams  of  energy  is  not  feasible  for  system  gecmetrles  applicable  to  this 
method . 


IV 


PROGRAM  FOR  MEXT  PHASE 


The  objective  of  the  next  phase  Is  to  design  experiments  capable  of 
verifying  the  assumptions  and  results  of  the  analyses  perfomed  in  this  phase ^ 
and  to  consider  the  instrumentation  required  to  conduct  the  experiments  de¬ 
signed.  Experiments  will  be  designed  to  determine  the  parameters  listed  be¬ 
low. 


A.  Scattering  From  Natural  Turbulence 

1.  'Rie  motion  of  turbulence  relative  to  the  wind. 

2.  The  dimensional  scale  of  low-level  turbulence. 

3.  fHie  refractive  index  variations  in  low-level  turbulence. 

4.  The  random  turbulent  velocities  and  their  contribution  to 
Doppler  broadening. 

5.  The  probability  of  occurrence  of  detectable  turbulence. 

B.  Electromagnetic  Scattering  From  Acoustic  Waves 

1.  Maximum  range  as  a  function  of  system  frequencies. 

2.  Reflection  coefficient  and  scattering  cross  section  of  acoustic 
waves  for  verification  of  theory  and  extension  of  analysis. 

3.  Attenuation  due  to  turbulent  scattering  of  the  acoustic  wave. 

4 .  Acoustic  beam  displacement  by  wind  and  turbulence . 

5.  Radar  scattering  from  turbulence  of  the  same  scale  as  the 
acoustic  wavelength. 

C.  Infrared  Tracking  of  an  Artificially  Heated  Bubble 

No  experimental  design  is  reconnended  for  this  method. 


-  27  - 


V.  PERSOWIEL 


The  following  personnel  have  contributed  to  the  project  for  the  num¬ 
ber  of  man-hours  listed. 

Harold  L.  Stout,  Assistant  Director,  Engineering  Division  -  32  man-hours. 

Paul  C.  Constant  I  Jr.,  Head,  Electronics  Engineering  Section  -  48  man-hovirs. 
Richard  W.  Fetter,  Senior  Engineer  -  400  man-hours. 

Paul  L.  Smith,  Jr.,  Senior  Engineer  -  222  man-hours. 

George  E.  Chambers.  Assistant  Engineer  -  32  man-hours. 

B«  L.  Jones,  Associate  Engineer,  100  man-hours,  B.S.  In  E.E.,  1947, 
University  of  Kansas,  Graduate  Studies  In  Physics,  1947-1948,  University  of 
Kansas.  Mr.  Jones'  special  fields  are  transmission  lines  and  electromagnetic 
radiation  and  propagation.  His  technical  experience  Includes  work  on  vibration 
suppression,  automation,  and  propagation  phenomena  with  regard  to  random  and 
Induced  disturbances  In  the  atmosphere. 

Harold  F.  Schick,  Senior  Engineer,  288  mun-hours,  B.S.  In  Electrical 
Engineering,  1954,  Cooper  Union;  M.S.  In  Electrical  Engineering,  1961,  Univer¬ 
sity  of  Arizona.  Mr.  Schick  has  thirteen  years' e]q>erlence  In  the  design  and 
analysis  of  electronic  systems.  The  past  seven  years  he  has  been  largely  con¬ 
cerned  with  air-borne  systems  having  stringent  requirements  with  regard  to  re¬ 
liability,  bandwidth  utilization,  minimum  oolse,  and  high  energy  density  per¬ 
formance.  A  significant  phase  of  the  design  and  analysis  has  entailed  evalua¬ 
tion  and  optimization  techniques. 

Robert  M.  Stewart,  Jr.,  Consultant,  56  man-hours,  B.S.  in  Electrical 
Engineering,  1945,  and  Ph.D.  In  Physics,  1954,  Iowa  State  College.  Mr.  Stewarth 
special  fields  are  thermal  and  mechanical  atmospheric  turbulence  and  digital 
and  cmalog  computer  solutions  of  oieteorologlcal  problems.  His  technical  ex¬ 
perience  Includes  teaching  electromagnetic  wave  theory. 


28  - 


BlBLIOQRAPHy 


1.  Papa,  R.  J.,  "Atmospheric  Turbulence,"  Oeophyslcs  Corp.  of  America,  Scien¬ 

tific  Report  No.  7,  for  Air  Force  Cambridge  Research  Laboratories  AD  257- 
642,  April,  1961. 

2.  Richards,  E.  J.,  and  J.  E.  Ff.  Williams,  "Some  Recent  Developments  in  Jet 

Noise  Research,"  University  of  Southampton,  Department  of  Aeronautics 
and  Astronautics,  Report  No.  118,  August,  1959.  See  edso  Bvill,  M.  K., 
"Instrumentation  for  and  PKlimlnary  Measurements  of  Space  -  Time  Correla¬ 
tions  and  Convection  Velocities..."  University  of  Southampton,  Depart¬ 
ment  of  Aeronautics  and  Astronautics,  Report  No.  149,  August,  1960. 

3.  Bricourt,  Pierre,  et  al.,  "Ultrasonic  Technique  for  the  Determination  of 

Wind  Velocity  Measurements/*  Emerson  Research  Laboratories,  1140  East- 
West  Highway,  Silver  Springs,  Maryland,  Final  Report  No.  FR-5122,  Con¬ 
tract  No.  NObs -64140  to  Bureau  of  Ships,  Department  of  the  Navy;  3  April 
1961. 

4.  Tillman,  R.  Q.,  et  ed.,  "A  Study  of  Clear-Air  Angels  by  the  Use  of  Hori¬ 

zontal  and  Vertical  Trajectories/'  Proc.  Ninth  Weather  Radar  Conference, 
Kansas  City,  Missouri,  pp.  241-246,  23-26  October  1961. 

5.  Stewart,  P.  0.,  "Radar  in  Turbulence  Research,"  ibid.,  pp.  206-211. 

6.  First  Phase  Report,  p.  44. 

7.  Smith,  P.  L.,  Jr.,  and  R.  W.  Fetter,  "Remote  Measurement  of  Wind  Velocity 

by  the  Electromagnetic  -  Acoustic  Probe.  I:  System  Analysis.  II;  Experi¬ 
mental  System,"  Conference  Proceedings,  Fifth  Natlonail  Convention  on 
Military  Electronics,  Washington,  D.  C.,  pp.  46-59,  June,  1961. 

6.  Fetter,  R.  W. ,  P.  L.  Smith,  Jr.  and  V.  Klein,  "Design  Concept  Study  for  a 
Wind  Profile  Instrument,"  Final  Report,  Contract  No.  DA-23-072-506-0RD-11 
to  U.  S.  Army  Rocket  and  Guided  Missile  Agency  from  Midwest  Research  In¬ 
stitute,  November,  1960. 

9.  Private  communication  from  Mr.  R.  W.  Sutton,  The  Boeing  Co.,  Seattle,  Wash. 

10.  Fetter,  R.  W.  and  P.  L.  Smith,  Jr.,  "A  Study  of  Possible  Applications  of 
the  EMAC  Probe  Aboard  Navy  Ships,"  Final  Report,  Contract  No.  NOw  61- 
0720-d  to  Bureau  of  Weapons,  Department  of  the  Navy  from  Midwest  Reseeurch 
Institute,  pp.  14-20,  October,  1961. 


-  29 


BmiOQRAPHY  (Concluded) 


11.  Stevens,  D.  W.,  "Cedar  Hill  Meteorological  Research  Facility,"  Instrumen¬ 
tation  for  Geophysics  and  Astrophysics,  No.  16,  GRD,  Air  Force  Cambridge 
Research  Laboratories,  ASTIA  No.  AD  263  847,  June,  1961. 


-  30  - 


APEEWDIX  A 


ANALYSIS  OF  EI£CTROMAGNETIC -ACOUSTIC  (EMAC)  PROBE 


This  analytical  treatment  of  the  EMAC  probe  was  originally  per* 
forme'’  "-"^tract  Ho.  DA-23-072 -506-0RD-11,  with  the  Army  Rocket  and 
Gul  -d  Missile  Agc...^  is  included  here  to  make  available  the  mathe¬ 

matical  procedures  and  resuxw-  xs  they  apply  to  the  present  analysis. 


I.  RAN(ai  ffiRFORMAMCE  OF  THE  EMAC  SYSTEM 


The  range  performance  of  the  EMAC  probe  for  remote  measurement  of 
wind  velocities  is  limited  chiefly  by  the  decrease  in  the  reflection  coef¬ 
ficient  of  the  acoustic  disturbance  with  increasing  distance  from  the  source. 
This  decrease  is  the  result  of  two  factors:  the  divergence  of  the  acoustic 
wave  and  the  absorption  and  scattering  by  the  atmosphere  of  some  of  the 
acovLstic  energy.  Since  the  absorption  Increases  with  the  acoustic  frequency^ 
the  use  of  as  low  a  frequency  as  is  compatible  with  other  reqid.rements  of 
the  system  is  indicated. 


A.  CalcxUation  of  Received  Power 


To  obtain  a  quantitative  estimate  of  the  ultimate  range  of  the 
EMAC  system,  consider  a  system  in  which  the  acoustic  and  radar  sources  are 
located  at  the  same  point.  This  arrangement  is  advantageous  because  the 
spherical  acoustic  wavefronts  act  like  concave  mirrors  in  focusing  the  re¬ 
flected  electromagnetic  energy  back  towsu^  the  source. 

1.  Reflected  power;  The  microwave  power  density  impinging  on 
the  acoustic  wave  train  at  a  distance  r  from  the  source  is 


Pi 


Pt  G 

4nr2 


(A-1) 


The  power  reflection  coefficient  p2  of  the  acoustic  wave  train  has  been 
c(»iputed  previously. V  The  reflected  power  per  unit  area  of  the  acoustic 
wavefront  will  be 


Pr 


Pi  X 


(A-2) 
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The  total  reflected  power  is  the  product  of  this  reflected  power  density  and 
the  area  of  the  acoustic  distuxhance  illuminated  by  the  radar  signal.  Be¬ 
cause  of  the  divergence  of  both  the  acoustic  auid  radar  beams,  the  illum¬ 
inated  area  will  Increase  in  proportion  to  the  square  of  the  distance.  The 


illuminated  eurea,  of  course,  never  can  be  greater  than 


and,  in 


general,  will  be  only  some  fraction  f  of  this  quantity.  Then  the  total 
reflected  power  will  be 


hi  p2  isft 

4TTr^  ® 


«  P^p2f  .  (A-3) 

The  value  of  f  depends  on  the  matching  of  the  geometry  of  the 
acoustic  and  radar  beams.  For  best  results,  the  beam  widths  should  be  made 
equal  —  this  will  Insure  that  most  of  the  radaur  power  is  concentrated  where 
the  acoustic  beam  is  most  intense,  hence,  where  the  reflection  coefficient 
is  greatest. 


2.  Received  power;  Because  of  the  focusing  action  of  the  acous¬ 
tic  wavefronts,  a  rather  large  fraction  g  of  the  total  reflected  power 
will  be  returned  to  the  radar  receiver.  The  total  received  power  will  be 


P,  «  P^pSfg  .  (A-4) 

The  product  fg  is  a  measure  of  beam  and  wavefront  alignment. 

The  deiJendence  of  the  power  reflection  coefficient  on  the  acoustic 
signal  may  be  expressed  asV 

p2  =  K  (hf  ,  (A-5) 

where  the  peurameter  K  has  the  value  61.6  for  an  acoustic  signal  10  cycles 
long  and  is  proportional  to  the  square  of  the  number  of  cycles  in  the  acous¬ 
tic  signal.  The  relationship  of  the  quantity  )  to  the  intensity  of 
the  acoustic  signal  Is^ 


-  1.2  X  10"^^  x\j^-  , 


Using  these  values,  the  reflection  coefficient  of  the  acoustic  signal  is 


p2  «  1.44K  X  10“®®  X 


(fe) 


(A-7) 


The  acoustic  intensity  decreases  with  distance  from  the  source  ac< 
cording  to  the  equation 


I  .  I 


i(t)“ 


(A-8) 


The  power  reflection  coefficient  is  therefore 

L.44k(^)' 


p2  «  1. 


2  „  I, 

X  10'2®  X  X  10 
■•■o 


{A-9) 


For  this  value  of  reflection  coefficient,  the  power  received  at  the  antenna 

) 

ri®(l.44)K  X  10“®®  X  I,  X  10  ^ 

Pr  -  Pt  - = - - -  f  6 


r"  Ir 


or 


p  r.  K  X  I.  X  10 

.Pfi  A  1 

~  =  1,44  X  10*^  - 


rP  ^o 


fg  . 


(A-10) 


B.  Typical  Values  of  Ttansmission  loss 

1.  Values  of  the  pturameters;  An  optimistic  estimate  of  the  re¬ 
ceived  power  is  obtained  from  the  following  values: 

(fg)  a  1;  =  140  decibels  at  r^  «  10  ft. 
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(this  value  of  Ii  has  been  attained  by  an  experimental  siren).  The  value 
of  K  depends  upon  the  number  of  cycles  employed  in  the  acoustic  signal 
which,  in  turn,  is  limited  by  the  remge  resolution  requirement.  For  pur¬ 
poses  of  the  subsequent  calculation,  an  acoustic  wave  train  length  of  25  ft. 
will  be  used.  It  is  believed  that  this  length  represents  the  maximum  allow¬ 
able  range  uncertainty,  \AiiIe  still  producing  usable  reflections.  At  a  fre¬ 
quency  of  22  kc.  (wavelength  1.5  cm.)  a  25-ft.  long  acoustic  signal  will 
consist  of  500  cycles  In  the  direction  or  propagation.  The  corresponding 
numbers  of  cycles  In  25  ft.  for  other  frequencies  are  Indicated  in  Table  A-I. 


TABLE  A-I 

VALUES  OF  THE  PARAMETER  K  FOR  SEVERAL  ACOUSTIC  FREQUENCIES 


Acoustic 

Frequency 


Acoustic 
Wavelength 
(cm. ) 


No.  of  Cycles 
in  25  Ft. 


K 


22  1.5 

11  3 

7.33  4.5 

5.5  6 


500 

154,000 

250 

38,500 

166 

17,000 

125 

10,000 

The  value  of  the  absorption  constant  a  depends  in  a  rather  com¬ 
plicated  way  upon  frequency  and  relative  humidity.^  Values  for  the  absorp¬ 
tion  constant  at  30,  50  and  70  per  cent  relative  humidity  are  plotted  as  a 
fxinction  of  frequency  in  Fig.  A-1;  the  data  were  taken  from  Knudsen's 
curves. 5/  At  50  per  cent  relative  humidity  the  variation  of  a  with  fre¬ 
quency  Is  approximately  linear,  with  the  equation 


a  =  (0.0107  fa  -  0.02)  db/ft 


(A-11) 


representing  the  values  for  frequencies  above  4  kc.  Note  that  for  other 
relative  humidities  the  dependence  of  a  on  freqviency  is  not  linear;  for 
lower  relative  humidities  the  absorption  will  be  higher  than  at  50  per  cent 
R.H.,  and  for  higher  humidities  it  will  be  lower. 

Knudsen's  data  extend  only  to  10  kc.  in  frequency;  the  values 
listed  below  for  the  absorption  constant  a  at  higher  frequencies  have 
been  obtained  by  extrapolating  the  linear  relationship  of  (A-11): 
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Fig.  A-1  -  Absorption  Constant  for  an  Acoustic  Wave  in  Air 
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I 

I 


a 

a 

a 

a 


0.22  db/ft  at 
0.10  db/ft  at 
0.06  db/ft  at 
0.04  db/ft  at 


22  ke,,  SO  per  cent  R.H. ; 
11  kc.,  50  per  cent  R.H.; 
7.3  kc.,  50  per  cent  R.H.; 
5.5  kc.;  50  per  cent  R.H. 


(A-12) 


A  plot  of  acoustic  Intensity  vs.  range  for  these  frequencies  is  shown  in 
Fig.  A-2. 


2.  Values  of  tranasdasion  loss;  Using  these  values  in  (A-10), 
the  following  equations  for  received  power  at  the  various  acoustic  frequen¬ 
cies  are  obtained: 


kc 


2.21  X  10’^  i2 


-0 


.022(r-10) 

"15 - 


kc 


5.55  X  10 


Q  ^Q.0.0l(r.l0) 

— ;5 — 


t 


kc 


2.45  X  10 


g  ^jj-0.006(r.l0) 
^2 


1.44  X  10 


.6  io'^-004(r~10) 


kc 


(A-13) 


The  above  ratios,  in  decibels  of  transmission  loss  where 


db  -  -10  log  —  , 


(A-14) 


are  plotted  as  a  function  of  range  in  Fig.  A- 3.  Assuming  a  rather  ideedlzed  re¬ 
ceiver  sensitivity  of  -135  dtan,  values  of  the  maximum  usable  transmission 
loss  for  transmitted  powers  of  10  watts,  100  watts  and  1  kw.  have  been  indi¬ 
cated  by  the  horizontal  lines  on  the  figure.  The  meaning  of  these  curves 
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DISTANCE  FROM  SOURCE  IK  FEET 


Pig.  A-3  -  Round  Trip  Transmission  Loss  vs.  Distance  for  Ideal  System 


will  be  clearer  from  the  following  example:  suppose  a  transmitter  power  of 
100  watts  is  used  with  a  receiver  of  sensitivity  -  135  dbm.  Then  for  various 
choices  of  acoustic  frequency  (and  the  corresponding  radar  frequency  deter¬ 
mined  from  Xq  >  2Xg^)  the  mayimm  range  will  be  as  given  in  Table  A-II. 


TABIE  A-II 


MAXIMUM  RAWffl!  OF  EMAC  SYSTEM  FOR  VARIOUS  FREftUEHCIES 


Acoustic 

Radar 

Maximum 

Frequency 

Frequency 

Range 

(kc.)  (nc»)  (ft.) 


22 

10,000 

400 

11 

5,000 

760 

7.3 

3,300 

1,140 

5.5 

2,500 

1,570 

C.  Maximum  Receiver  Sensitivity 

The  available  noise  power  at  the  input  to  the  receiver  will  be  of 
the  order  of  4  x  10“^  watte  per  megacycle  of  receiver  bandwidth. If  a 
receiver  with  a  noise  figure  of  3  db.  is  used,  the  equivalent  noise  power 
input  to  the  receiver  will  be  6  x  10*^  watts  per  megacycle  of  bandwidth. 

The  required  bandwidth  depends  on  the  amount  of  Doppler  shift  vAklch  will  be 
contributed  by  the  wind  velocity;  for  maximum  winds  of  75  mph,  the  required 
bandwidth  will  be  roughly  150  times  the  Doppler  shift  in  cycles  per  second 
per  mile  per  hour  of  wind  velocity.  Typical  values  of  receiver  bandwidth, 
equivalent  noise  power  input  and  maximum  attainable  receiver  sensitivity  for 
a  receiver  capable  of  detecting  a  signal  equal  in  power  to  its  equivalent 
input  noise  power,  are  shown  in  Table  A-III. 


D.  Scattering  Due  to  Tuibulence  in  the  Atmosphere 

An  additional  factor  which  may  affect  the  probe  performance  at 
longer  ranges  is  the  scattering  of  the  acoustic  signal  by  atmospheric  turbu¬ 
lence.  This  scattering  will  be  seen  by  the  system  as  Increased  attenuation 
of  the  acoustic  wave.^  The  amount  of  scattering  attenuation  depends  upon 
both  the  scale  of  the  turbulent  eddies  and  the  magnitude  of  the  fluctuations 
in  the  properties  of  the  atmosphere.  Because  of  the  lack  of  Information  on 
these  effects,  it  is  difficult  to  establish  a  quantitative  estimate  of  the 
attenuation  constant  to  be  expected. 
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TAB1£  A-III 

MAXIMUM  AITAINABIE  8EMSITIVITY  FOR  BECEIVER  WITH  5  DB.  H0I3E  FIOURE 


Acoustic 

Frequency 

(kc.) 

Microwave 
Frequency 
(me. ) 

Doppler 

^ift 

(cps/mph) 

Approximate 

Bandwidth 

Required 

(kc.) 

Equivalent 

Noise  Power 
Input 

(x  lO"!"^  watts) 

Maximum 

Receiver 

Sensitivity 

(dbm) 

22 

10,000 

30 

55 

4 

-134 

11 

5,000 

15 

2.5 

2 

-137 

7.3 

3,300 

10 

2 

1.6 

-138 

5.5 

2,500 

7.5 

1.5 

1.2 

-139 

However,  analysis  shows^  that  It  Is  probable  that  the  attenuation  of  the 
acoustic  signal  due  to  scattering  could  exceed  that  due  to  atmospheric  ab> 
sorption.  The  range  performance  of  the  EMAC  prbbe  would  be  reduced  by  this 
additional  attenuation,  but  no  convenient  quantitative  evaluation  of  the  re¬ 
duction  is  possible. 


The  attenuation  constant  due  to  scattering  is,  however,  known  to  be 
proportional  to  the  square  of  the  acoustic  frequency.^  The  use  of  a  lower 
acoustic  frequency  will  therefore  reduce  the  attenuation  from  this  cause  as 
well  as  frcm  absorption. 


E.  Experimental  Verification 


Experimental  measurements  were  made  on  the  basic  EMAC  probe  to 
verify  the  theoretically  derived  equation  for  the  function,  Pj./Pt  t  for  the 
system,  based  upon  an  assumed  value  for  the  factor  fg  in  (A-10).  For 
reference,  the  basic  equation  is: 


Pr 

Pt 


1.44  x  10 


•26  ‘1 


T 

r,  K  X  I,  X  10 


fG  , 


r2  I. 


(A-10) 


where  the  parameters  and  their  measured  values  are: 

Pp  =  power  received  at  antenna  =  -96  dbm  (receiver  sensitivity); 
Pt  =  power  transmitted  =  +34.0  dbia  (3  watts); 
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Iq  ■  reference  acoustic  intensity  (defined)  10"^®  watts/cm^  =  0  db.; 
>  acoustic  intensity  at  a  range  r]^  »  140  db.; 

-  140  db.  or  10^^; 

lo 

r^  «  range  at  \diich  is  measured  -  10  ft.; 

r  «  maximum  range  of  detection  ■  90  ft.;  and 

a  -  the  acoustic  absorption  constant  «  0.22  db/ft  for  22  kc. 

The  parameters  \dilch  %fere  not  measured  are: 

f  «  the  fraction  of  the  acoustically  disturbed  area  which  is 
lllvnolnated  by  microwaves; 

g  «  the  fraction  of  the  reflected  energy  which  is  intercepted  by 
the  receiving  antenna;  and 

K  ■  derived  constant  determined  by  length  of  acoustic  wave  train 
«  10^  for  1^  wavelengths  in  a  e-ms.  wave  train. 

The  measurement  of  f,  g,  and  K  for  the  experimented  system  was 
beycmd  the  scope  of  this  program,  so  complete  verification  of  (A-IO)  was  not 
possible.  To  complete  the  calculation,  the  derived  value  of  K  *  10^  etnd 
an  Initial  assumption  of  fg  >  1  were  used.  The  assvmptlon  for  fg  implies 
that  the  system  is  ideal,  that  is,  all  microwave  energy  interacts  with  tne 
acoustic  disturbance  and  all  energy  that  is  reflected  is  focused  on  th.-  re¬ 
ceiving  antenna.  This  optimistic  assumption  gives  a  predicted  maximum  range 
of  142  ft.  VThen  compared  with  the  measured  value  of  90  ft.,  a  discrepancy  of 
16  db.  is  indicated  in  the  over-all  transmission  loss.  This  is  shown  in 
Fig.  A-4,  \diich  is  a  plot  of  Pj/P^  in  decibels  as  a  function  of  range  for 
the  experimental  system.  A  more  realistic  estistate  of  fg  was  obtained  by 
assuming  that  all  the  reflected  energy  was  returned  frcm  interaction  within 
the  5*,  half -power  beam  width  of  the  main  lobe  of  each  antenna  and  that  the 
acoustic  wave  front  geometry  fell  half  way  between  the  siberlcal  shape  pro¬ 
duced  by  a  point  source  and  the  plane  waves  produced  by  a  perfectly  fo.used 
reflector.  From  these  considerations,  a  value  of  fg  «  l/30  was  obtained. 
While  this  number  is  not  exact,  it  is  much  more  realistic  than  the  first  as¬ 
sumption  and  leads  to  a  predicted  maximum  range  of  93  ft.,  or  a  discrepancy 
of  1.5  db.  in  transmission  loss.  This  agreement  is  sufficiently  good  to 
Justify  using  the  theory  to  predict  the  performance  of  systems  operating  at 
other  frequencies  and  power  levels. 
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Fig.  A-4  -  Round  Trip  Transmission  Loss  vs.  Distance  for  Experimental  System 


To  obtain  complete  correlation  between  the  theoretical  and  experl> 
mental  system  characteristics,  accurate  measurements  of  f  ,  g  ,  and  K 
should  be  made,  preferably  at  lower  system  frequencies  where  much  greater 
ranges  are  predicted  and  the  system  performance  should  be  more  sensitive  to 
parameter  variations. 


^VSTOM  COMFIQURATIOilS  FOR  QBTMMING 
-"^^noiTY  PROFimS 


A.  Requirements  of  the  S^tem 

The  EMAC  probe  meaaxires  the  vector  svr:  of  the  velocity  of  propaga¬ 
tion  of  the  acoustic  signal  through  air  eind  the  component  of  wind  velocity 
In  the  direction  of  propagation.  Since  the  wind  velocity  may  be  represented 
as  a  vector  with  three  spatial  components,  the  system  will  be  required  to 
gather  four  Independent  pieces  of  data  from  which  the  four  quantities  v^  , 
Wx  ,  Wy  and  W2  may  be  determined.  In  general,  these  four  pieces  of  data 
will  be  obtained  from  four  soundings  made  by  an  EMAC  probe  or  probes  In  dif¬ 
ferent  directions. 

In  particular  situatlcms  only  three  such  soundings  may  be  required. 
This  would  be  the  ease  If  only  horizontsil  components  are  desired,  or  if  It 
can  be  assumed  that  the  vertical  component  of  the  wind  is  much  smaller  them 
the  other  two  components.  In  the  usuad  meteorological  reports  the  vertical 
component  Is  not  even  considered.  If  the  vertical  component  Is  neglected, 
only  three  soundings  are  required,  since  only  three  of  the  unknowns  remain. 


B.  System  Using  One  Steerable  EMAC  Probe 

For  the  case  of  low  elevation  angles,  the  two  horizontal  components 
of  the  wind  velocity  vector  can  be  obtained  by  taking  three  soundings  in  dif¬ 
ferent  directions  with  a  single  EMAC  probe  operating  from  a  fixed  site  at  the 
launching  point.  The  directions  of  the  soundings  are  shown  schematically  In 
Fig.  A-5. 


tor  is 


1.  Method  of  obtaining  wind  velocity  components:  If  the  wind  vec- 


W  *  Wx  1  +  Wy  J  +  Wz  k  , 


(A-15) 
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Fig.  A-5  -  Sounding  Geometry  of  Single  Steerable -Probe  System 


then  the  soundings  give  velocity  Information  which  may  be  represented  by 

a*VB+W*A  ; 

b  «  Vs  +  W  •  B  ;  and 

c  «  Vs  +  W  .  C  .  (A-16) 

The  unit  vectors  are  obtained  by  referring  to  the  diagram  of  Fig.  A-5. 

A  B  -cos  0  sin  6  1  cos  0  cos  6  J  sin  k  ; 

B  B  cos  J  sin  0  k;;  and 

C  =  cos  0  sin  01+  cos  0  cos  0  J  +  sin  0  k  .  (A-17) 

Substituting  these  values  Into  (k-16), 

a  B  Vg  -  Wjj  cos  0  sin  B  +  Wy  cos  0  cos  0  +  sin  0  ; 

b  =  Vg  +  Wy  cos  0  +  Wj  sin  0  ;  and 

c  B  Vg  +  Wjj  cos  0  sin  ©  +  Wy  cos  0  cos  0  +  sin  0  .  (A-18) 

Since  U2  Is  small  and  the  sines  of  the  elevation  angles  are  also 
smsill,  the  W2  terms  may  be  neglected,  leaving  the  result: 

a  =  Vg  -  Wjj  cos  0  sin  0  +  Wy  cos  0  cos  0  ; 

b  =  Vg  +  Wy  cos  0  ;  and 

c  =  Vg  +  Wx  cos  0  sin  +  Wy  cos  0  cos  0  .  (A-19) 

This  set  of  three  simultaneous  equations  may  be  solved  to  give  the  wind 
velocity  components  as  follows: 


f  =  c  -  a 
^  2  cos  0  sin  0 

I  _  a  +  c  -  2b _ 

^  2  cos  0  (cos  0-1) 


(A-20) 
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2.  Errors  In  the  calculated  wind  velocity  ccmponents;  In  e*»ieral, 
each  sounding  results  In  a  reading  (a,  b  or  c)  uhlch  Is  In  error  by  some 
amount  e  .  The  elevation  angle  0  is  small,  so  its  cosine  is  nearly 
equal  to  unity.  The  errors*  in  the  calculated  wind  velocity  components  axe 
then: 


A(Wx)  a 

1  ^  1 

|sin  91 

A(Wy)  a 

2e 

cos  9 

(A-21) 


The  values  of  these  errors  are  plotted  as  functions  of  0  in  Fig.  A- 6.  The 
error  in  the  range  component  Wy  is  considerably  larger  than  that  in  the 
cross-wind  component  Wx  .  The  significance  of  the  curves  may  be  Illustrated 
by  an  example:  suppose  A  »  45*  and  e  «  0.5  mph.  Then  the  errors  in  the 
respective  wind  velocity  components  are 


A(Wx)  <  1.41  e  (approx.  0.7  mph)  ;  and 


A(Wy)  <  6.65  e  (approx.  3.4  mph)  . 


Soundings  with  an  accureu:y  of  1  0.5  mph  require  determining  the  Doppler  shift 
to  about  1  part  in  1,500,  since  the  velocity  of  propagation  of  the  sound  wave 
contributes  about  750  mph  to  the  velocity  being  sensed. 

An  additional  error  is  contributed  by  the  terns  in  which  were 
neglected  in  solving  (A-18).  This  is  more  than  Just  a  question  of  an  unknown 
vertical  wind  component,  for  the  presence  of  the  vertical  coaponent  intro¬ 
duces  an  error  into  the  calculated  values  of  the  other  components  if  (A-20) 
is  used.  Studlesi/  have  shown  that  the  vertical  component  of  wind  velocity 
may  in  some  cases  be  as  hi£^  as  5  mph.  In  such  a  situation,  the  neglected 
terms  have  values  of  about  1  mph  for  an  elevation  angle  of  12*.  This  will 
set  a  lower  limit  on  the  error  occurring  in  these  measurements,  unless  a 
fourth  sounding  is  used  so  that  may  also  be  calculated.  Unfortunately, 
the  likelihood  of  encountering  such  verticeJ.  wind  velocities  at  low  altitudes 
is  not  known  accurately. 


*  If  the  errors  in  the  soundings  were  random,  the  errors  in  the  calculated 
wind  components  could  be  determined  as  the  square  root  of  the  sum  of  the 
squares  of  the  errors  in  the  individual,  soundings.  At  present,  however, 
there  is  no  basis  for  expecting  random  errors,  so  simple  addition  of 
errors  has  been  used. 
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Fig.  A-6  -  Errors  in  Wind  Velocity  Components  Determined  from  Three 

Soundings  with  a  Single  EMAC  Probe 
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In  the  above  anaXyels,  data  from  the  three  soundings  were  combined 
so  as  to  eliminate  the  Vg  terms  of  (A>19).  Inspection  of  (A«19)  for  b  , 
however,  shows  that  If  a  good  estimate  of  Vg  Is  available,  the  error  In 
Wy  would  be  equal  to  the  sum  of  e  and  the  error  In  Vg  .  Under  such  a 
condition,  error  In  range  component  Is  the  same  as  that  discussed  later  In 
Section  II>D-2  of  this  Appendix. 

3.  Discussion  of  the  slngle»prdbe  system;  The  single -probe  sys¬ 
tem  has  the  obvious  advantage  of  requiring  only  one  EMAC  probe.  In  contrast 
to  the  mviltl-probe  systems  described  In  Sections  II-C,  D  and  E.  In  addi¬ 
tion,  no  provisions  for  scanning  the  flight  path  with  the  single  probe  are 
required.  Each  of  the  three  soundings  Is  made  with  a  fixed  azimuth  angle 
(determined  by  the  choice  of  9  )  and  a  fixed  elevation  angle  (dependent 
upon  the  flight  path  chosen). 

The  major  objection  to  the  use  of  the  single -probe  system  arises 
because  (as  can  be  seen  from  Fig.  A-6)  large  values  of  9  must  be  used  If 
good  accuracy  Is  to  be  obtained.  Consequently,  the  data  gathered  by  sound¬ 
ings  A  and  C  come  from  points  rather  far  removed  from  the  desired  loca¬ 
tion.  For  example.  If  9  «  30*  Is  selected  as  giving  sufficient  accuracy  In 
the  calculated  wind  velocity  components;  then  for  a  point  at  range  r  from 
the  probe,  the  wind  velocity  components  will  be  calculated  from  data  ob¬ 
tained  at  one  point  on  the  path  and  two  other  points  at  a  distance  1  ,  from 
the  first  one,  where 


1  »  2r  sin  15*  ,  (A-22) 

as  Illustrated  in  Fig.  A-7.  When  r  >  1,000  ft.,  then  d  *  520  ft.;  so  that 
two  of  the  soundings  are  made  at  points  separated  by  1,000  ft.  While  (A-18) 
contains  the  Implicit  assumptions  that  both  the  wind  velocity,  W  ,  and 
acoustic  propeigatlon  velocity,  Vg  ,  remain  constant  during  the  time  Interval 
In  which  the  soundings  are  made  and  have  the  same  value  for  all  three  sound¬ 
ing  points,  it  appears  quite  likely  that  with  separations  as  great  as  the 
range,  r  ,  these  assumptions  may  not  be  valid.  A  statistical  approach  might 
give  useful  Information  over  a  given  time  lntez*val  but  has  not  been  considered 
in  this  study.  The  distance  between  outer  soundings  could  be  reduced  by  de¬ 
creasing  the  angle  9  ,  but  the  errors  in  calculated  wind  velocities  will 
Increase  as  shown  In  Fig.  A-6. 

If  the  relatively  low  accuracy  of  the  single-probe  system  can  be 
tolerated,  the  system  would  be  attractive  because  of  its  simplicity.  The 
multi-probe  systems  described  below  require  rather  csureful  siting  which  may 
be  difficult  to  accomplish  under  many  conditions. 


48 


FLIGHT  HATH 


C.  System  Using  One  Fixed  and  Tvo  Steerable  EMAC  Probes 

A  second  method  of  obtaining  the  tvo  horizontal  wind  velocity  com> 
ponents  consists  of  taking  the  three  soundings  with  three  different  EMAC 
probes  operating  from  Independent  sites.  A  configuration  for  such  a  system 
is  shown  schematically  in  Fig.  A-6;  the  location  of  the  two  side  probes  B 
and  C  at  a  distance  r/2  down  the  path  represents  an  obvious  but  rather 
arbitrary  choice. 

1.  Method  for  obtaining  wind  velocity  components;  For  measure¬ 
ments  made  when  the  wind  vector  at  some  point  is  7  >  (A-15)^  the  soundings 
give  velocity  infoxnation  which  is  represented  by  (A-16).  In  the  three- 
probe  case>  the  unit  vectors  may  be  obtained  with  reference  to  Fig.  A-9  as: 


A  B  cos  ^  i  *  sin  ^  k  ; 

B  «  cos  Y  cos  ?  i  +  cos  Y  sin  9  J  +  sin  Y  k  ;  and 

C  ■  -  cos  Y  cos  01+  cos  Y  sin  0  J  +  sin  Y  k  .  (A-23) 

Substituting  these  values  into  (A-16), 
a  »  Vg  +  Wy  cos  0  +  Wg  sin  0  ; 

b  ■  Vg  +  Wjj  cos  Y  cos  P  +  Wy  cos  y  sin  0  +  sin  y  >  and 

c  “  Vg  -  Wjj  cos  Y  cos  0  +  Wy  cos  y  sin  0  +  W^  sin  y  .  (A-24) 

Again,  for  small  values  of  W^  and  small  angles  of  elevation,  the 
W2  terms  may  be  neglected.  This  leaves 

a  =  Vs  +  Wy  cos  0  ; 

b  =  Vg  +  Wx  cos  Y  cos  P  +  Wy  cos  y  sin  0  ;  and 
c  =  Vs  -  Wx  cos  Y  cos  0  +  Wy  cos  y  sin  0  .  (A-25) 

These  equations  may  be  solved  to  give  the  wind  velocity  components  as  follows; 
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Wg,  A-9 


Coordinate  System 


for  Ihree -Probe 


Configuration 
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w. 


b  -  c 


2  cos  Y  cos  0 


W 


y 


2(eos  Y  sin  9  -  cos  0) 


(A-26) 


2.  Errors  In  the  calculated  wind  velocity  components;  Each  sound 
ing  yields  a  reading  (a,  b  or  c)  which  is  in  error  by  some  amount  e  .  The 
cosines  of  the  elevation  angles  ,  y)  are  nearly  equal  to  unity  If  the 
angles  are  small.  Therefore,  the  exrors  in  the  calculated  wind  velocity  com 
ponents  are  approximately: 


A(Wx)  s 

e 

cos  6 

A(Wy)  £ 

2e 

sin  9 

(A-27) 


The  values  of  these  errors  sire  plotted  as  functions  of  0  in  Fig.  A>10;  In 
this  system  the  value  of  9  is  not  fixed,  but  must  be  swept  through  the 
lange  -45*  to  <»45*  to  cover  the  entire  flight  path.  This  must  be  accca* 
plished  by  scanning  the  flight  path  with  preibes  B  and  C  . 


The  Intepretation  of  these  curves  is  Illustrated  by  the  following 
example.  Assume  that  a  low-angle  path  is  used,  and  that  e  «  0.5  mph  ;  the 
errors  in  the  calculated  wind  velocity  components  at  various  points  of  the 
path  are  given  in  Table  A-IV. 


Position 
Launch  Point 
Hid -Range,  £ 
End  of  Rcmge, 


TABI£  A-IV 

TYPICAL  ERRORS  IH  WIKD  VELOCITY  COMPOHEMTS 
DETERMIHED  BY  A  THREE -IHOBE  SYSTEM 

(9  *  -45*)  1.41  e  0.7  mph 

(O  =  0*)  1.0  e  ft;  0.5  mph 

r  (0  =  +45*)  1.41  e  ft.-  0.7  mph 


A(Wy) 


1.17  e  ft:  0.6  mph 
2  e  ft:  1.0  mph 
6.85  e  ft:  3.4  mph 
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EMM  IN  WIND  VELOCITY  COIffOMT  (IN  WITl  OF  •) 


AZIMUTH  ANOLE  0 


Fig.  A-10  -  Errors  in  Wind  Velocity  Components  Deteindned  by 

Three -Probe  Sycten 
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As  in  the  case  of  the  slngle-prcibe  system^  an  error  Is  contrlbvtted 
when  the  terms  In  W2  in  (A>24)  are  neglected.  The  same  considerations  as 
to  the  magnitude  of  this  error  apply  in  the  three 'prcbe  case. 


3.  Relationship  between  sounding  angles;  The  elevation  angle  y 
will  be  dependent  upon  the  other  two  angles  9  and  ^  .  With  the  aid  of 
Fig.  A-9  the  relationship  among  these  angles  may  be  found  to  be 


Y  “ 


(sin  e  cos 


(A-28) 


Provision  must  be  made  when  scanning  the  flight  path  with  prdbes  B  and  C 
to  maintain  this  relationship  between  the  elevation  angle  y  and  the  azimuth 
angle  9  . 


4.  Discussion  of  the  three-probe  system;  The  arrangement  using 
three  EMAC  probes  is  attractive  for  the  measurement  of  wind  profiles  along 
low-elevation  trajectories  because  all  the  data  used  in  calculating  the  wind 
velocity  components  at  a  point  come  from  measurements  at  that  point.  No 
assumptions  concerning  uniformity  of  the  wind  velocity  over  extended  dis¬ 
tances  are  necessary,  in  contrast  to  the  situation  in  the  case  of  the  single - 
probe  system. 


D.  Substitution  of  Direct  Temperature  Measurements  for  One  of  the  Soundings 

If  the  velocity  of  sound  is  assumed  to  be  constant  over  the  path 
for  which  the  wind  velocity  profile  is  to  be  detemlned,  then  its  value  may 
be  determined  by  a  measurement  at  one  point  along  the  path.  Only  two  ad¬ 
ditional  pieces  of  data  are  needed  to  calculate  the  wind  velocity  components; 
these  may  be  obtained  by  soundings  made  with  only  two  EMAC  probes. 

The  only  variable  ^ich  has  a  significant  effect  on  the  velocity  of 
sound  in  air  under  usual  conditions  is  the  temperature  of  the  air. 2/  The 
rate  of  change  of  the  velocity  of  sound  with  respect  to  temperature  change  is 
about  0.75  mile  per  hour  per  degree  Fahrenheit.  Since  the  rate  of  change  of 
air  temperature  with  height  is  normally  about  0.55*F  per  100  ft.§/  for 
heights  not  too  close  to  the  ground,  the  error  resulting  from  the  assumption 
of  constant  temperature  along  the  path  will  not  preclude  its  application  for 
low-elevation  emgles  where,  at  the  maximum  range  of  the  EMAC  system,  the 
height  of  the  sounding  is  only  a  few  hundred  feet.  The  mesisurement  of  ed.r 
temperature  may  be  made  at  the  probe  site  by  a  thermocouple  or  other  suitable 
element,  or  the  local  sound  velocity  may  be  measured  directly  by  electronic 
means . 
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1.  Method  of  obtaining  wind  velocity  ccmponents;  To  obtain  the 
two  required  soundings,  only  two  EMAC  probes,  say  A  and  B  of  Fig.  A-8, 
need  be  used.  The  velocity  Information  obtained  by  these  two  probes  has  been 
previously  calculated  in  (A-25): 


a  =  Vs  +  Wy  cos  ^  ; 

b  a  Vs  +  Wx  cos  Y  cos  Q  +  Wy  cos  Y  sin  0  .  (A-29) 


The  third  "equation"  gives  Vg  directly  in  terns  of  the  measured  tempera¬ 
ture;!/ 


Vg  *  (717.3  +  0.75  T)  mph  . 


(A-30) 


The  wind  velocity  components  calculated  from  these  equations  are: 


cos  Y  cos  0 


tan  e  ^  ..  /tan© 

-  a - T  +  Vs  - 1 

V^cos  p 


cos 


cos 


i _ . 

Y  cos  ©y  ' 


and 


a  -  v« 


cos 


T  • 


(A-31) 


2.  Errors  in  the  two-probe  method;  In  addition  to  the  error  e 
in  the  soundings,  an  error  6  in  the  measured  velocity  of  sound  in  intro¬ 
duced  by  the  assumption  that  the  temperature  (and  hence  Vg)  is  constant 
along  the  flight  path.  If  the  cosines  of  the  elevation  angles  P  and  y 
are  approximately  equal  to  unity,  then  the  approximate  errors  in  the  ve¬ 
locity  components  are: 


e  (1  +  Isin  Pl)j 

j.  6(1  +  sin  ©) 

cos  0 

cos  P 

and 


i(Wy) 


s  e  +  6 


(A-32) 


Errors  in  the  temperature  sensing  instrument  are  included  in  6  . 
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The  size  of  the  errors  Involved  may  be  Illustrated  by  means  of  an 
example:  suppose  that  e  «  0.5  mph,  and  that  the  error  Introduced  as  6 
corresponds  to  2*F,  or  1.5  mph.  The  errors  In  the  calculated  wind  velocity 
components  along  the  path  are  shovm  In  Fig.  A-11.  The  errors  introduced  by 
neglecting  the  vertical  canponent  of  the  wind  velocity  have  not  been  In¬ 
cluded. 


3.  Discussion  of  the  tvo-probe  system;  The  two-probe  system  of¬ 
fers  a  modest  Improvement  in  simplicity  over  the  three -probe  system.  How¬ 
ever,  the  assumption  of  constant  temperature  along  the  flight  path  results 
In  larger  errors  In  the  calculated  wind-velocity  components.  Whether  the 
sacrifice  In  accuracy  offsets  the  reduced  complexity  of  the  system  must  be 
subjectively  decided. 

The  use  of  a  sound  velocity  measurement  to  replace  one  of  the 
soundings  may  be  particularly  attractive  In  adapting  the  three -probe  system 
to  tactical  applications. 
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APPENDIX  B 


ANALYSIS  OF  A  HEATED  BUBBIE  OF  GAS  IN  THE  FREE  ATMOSPHERE 


I.  THERMODYNAMIC  BEHAVIOR 


A  quantitative  assessment  of  the  behavior  of  a  bubble  of  heated  gas  In  the 
atmosphere  Is  undertaken  below  In  terns  of  the  simplest  form  of  modeling.  The 
bubble  Is  to  be  released  at  ground  level  and  Its  motion,  shape,  temperature 
and  density  are  determined  up  to  1500  meters.  A  simple  derivation  of  each  re¬ 
lationship  Is  presented  as  well  as  a  numerlcail  evaluation.  There  are  three 
treatments  for  each  relationship  corresponding  to  assumed  tenperature  lapse 
rates  for  the  atmosphere. 

As  stated  before,  the  most  simple  form  of  model  Is  used  for  this 
Investigation  In  order  to  provide  estimates  of  the  thermodynamic  quantities. 
Individual  assumptions  emd  simplifications  can  certainly  be  questioned  with 
respect  to  their  validity  In  the  real  world;  however  It  Is  felt  that  meaning¬ 
ful  conclusions  can  be  drawn  with  respect  to  feasibility  of  infrared  tracking 
of  the  heated  bubble  as  a  consequence  of  this  study. 


A.  Assumed  Model 

The  essential  description  of  the  model  is  given  below. 

1.  Turbulent  and  rooleculeur  diffusion  effects  on  the  shape  or  energy 
content  of  the  bubble  are  neglected. 

2.  Radlatlonal  exchange  between  the  bubble  and  Its  environment  are 
not  permitted. 


3.  The  thermodynamic  processes  occurring  in  the  bubble  are  strictly 
adiabatic  In  nature. 

4.  The  atmosphere  is  in  static  equilibrium  and  exerts  no  shear  or 
drag  forces  on  the  bubble  which  can  modify  its  shape  or  influence  its  motion. 

5.  The  atmosphere  can  be  characterized  by  a  uniform  teirperatiu:*e 
gradient.  In  this  investigation  there  are  three  gradients  assumed:  isothermal, 
-I'C/lOO  meters,  and  -2“C/l00  meters. 
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B.  E(iviatlon8  Used  In  Derivations  of  the  Ihemodynamlc  Relatlonshlpa 


iirtiere 


vrtiere 


where 


The  equations  that  will  be  used  In  this  derivation  are : 

1.  The  hydrostatic  equation  for  the  atmosphere: 

p  is  pressure, 

Z  is  height, 

g  is  acceleration  of  gravity,  and 
p  is  density  of  air. 

2.  The  perfect  gas  lav: 

p  =  RpT  ,  (B-2) 

p  is  pressure, 

R  is  universal  gas  constant/molecular  weight, 
p  is  density,  and 
T  is  absolute  temperature. 

3.  The  first  law  of  thermodynamics: 

dQ  =  Cy  dT  +  A  pd  ,  (B-3) 

dQ  is  amount  of  heat  exchanged  across  boundaries  of  system/unit  mass, 
Cy  is  specific  heat  at  constant  voliaue, 
dT  is  differential  temperature  change. 
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vrtiere 


A 

P 


is  the  reciprocal  of  mechanical  equivalent  of  heat, 
is  pressure,  and 

is  the  differential  change  in  volume. 
Characterization  of  the  free  atmosphere: 

®  Tao  “  cZ  , 

is  temperature  of  the  atmosphere  at  any  point  in  “K, 
is  temperature  of  the  atmosphere  at  Z  =  0  , 
is  the  lapse  rate  in  ‘C/lOO  meters,  and 
is  the  height  in  meters. 


(B-4) 


C.  Derivations  of  Thermodynamic  Relationships 

1.  Determination  of  bubble  temperature :  The  differential  of  (B-2) 
is 


pd  =  RdT  -  .  (B-5) 


Substituting  (B-5)  into  (B-3)  yields 


dQ  =  (Cy  +  AR)dT  -  APT  ^  .  (B-6) 

P 


Also,  AR  =  Cp-Cy  where  Cp  is  specific  heat  at  constant  pressure,  hence, 

dQ  =  Cp  dT  -  (Cp-Cy)  T  ^  .  (B-7) 

P 


Since  Y  =  -£  and  dQ  =  0  for  adiabatic  case,  (D-7)  becomes 
Cy 
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T  Y  P 


(B-0) 


Substituting  P  ==  Ip  f*”®®  (B-2)  in  (B-l)  yields 


^  = 


P 


(B-9) 


Substituting  (B-9)  in  (B-8)  and  assuming  P^  “  P^  • 


d_Tt  Vb-1 
Tb  "  ■  Vb 


(B-10) 


\rtiere  Tb  is  temperature  of  the  bubble. 

For  an  isothermal  atmosphere,  T  o  T  =  constant,  the  solution  of 

(B-10)  is 


For  finite  lapse  rate,  Tq  =  -  cZ  ,  the  solution  of  (B-10)  is 


Tb  =  T|jq 


(B-12) 


The  ratio,  Tb/Tbo  >  of  bubble  temperature  at  height  Z  to  its 
temperature  at  the  ground  is  tabulated  in  Table  B-I  for  various  heights  and 
lapse  rates,  a  . 
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TABLE  B-I 


RATIOS  OF  BUBBLE  TEMPERATURES 
Z(meter8) 


a("C/l00  meters) 

0 

200 

400 

600 

800 

1000 

1200 

1400 

1500 

0 

1 

0.999 

0.9S6 

0.981 

0.974 

0.967 

0.962 

0.954 

0.952 

1 

1 

0.993 

0.967 

0.98 

0.973 

0.967 

0.96 

0.953 

0.95 

2 

1 

0.994 

0.987 

0.98 

0.973 

0.966 

0.96 

0.953 

0.95 

2.  Determination  of  atmospheric  density:  The  differential  of  (B-2) 
Is 


dp  =  R(Tdp  pdT) 


(B-13) 


Substituting  (B-13)  into  (B-l)  yields 


(B-14) 


For  an  isothermal  atmosphere,  T  =  Toq  and  ^  =  0  ,  the  solution 


of  (B-14)  is 


az 


(B-15) 


For  finite  lapse  rate,  T^  =  T^o  -  r-Z  ,  the  solution  of  (B-14)  is 


f’a  “ 


(B-16) 


The  ratio,  Pa/Pao  >  atmospheric  density  at  height  Z  to  its 
density  at  the  ground  is  tabulated  in  Table  B-II  for  various  heights  and  lapse 
rates,  a  . 
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TABI£  B-II 


RATIOS  OF  ATMOSPHERIC  DENSITIES 

ZCmeters) 


a(*C/lOO  meters) 

~0~ 

200 

400 

600 

600 

1000 

1200 

1400 

1500 

0 

1 

0.978 

0.957 

0.935 

0.914 

0.893 

0.873 

0.853 

0.644 

1 

1 

0.963 

0.969 

0.953 

0.936 

0.923 

0.907 

0.691 

0.684 

2 

1 

0.991 

0.961 

0.972 

0.963 

0.953 

0.943 

0.934 

0.929 

3.  Determlnatlai  of  density  of  bubble:  It  is  assumed  that  at  all 
levels  .  Prom  (b-2), 

Pa  Ra  Ta  “  Pb  Tb  ; 
or 


For  an  Isothermal  atmosphere,  substitution  of  (B-ll)  eund  (B-15)  into 
(B-17)  yields 


^a  ^ao 

K 


(■ 


^ao^a^^ 


jl 


But 


Q  _  "^ao  0 


} 


therefore 
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(B-18) 


For  finite  lapse  rate,  substitution  of  (B-12)  and  (B-16)  into  (B-17) 

yields 


Pb  =* 


(B-19) 


The  ratio,  P^/Pbo  *  bubble  density  at  hei^t  Z  to  its 

density  at  the  ground  is  tabulated  below  in  Table  B>III  for  various  heights 
and  lapse  rates,  a  . 


TABLE  B-III 

RATIOS  OF  BUBBLE  DENSITIES 
ZCmeters) 


a(*C/lOO  meters) 

0 

200 

400 

600 

§00 

1000 

1200 

1400 

1500 

0 

1 

0.985 

0.  969 

0.954 

0.938 

0.925 

0.908 

0.895 

0.887 

1 

1 

0.983 

0.969 

0.953 

0.936 

0.923 

0.907 

0.891 

0.884 

2 

1 

0.984 

0.967 

0.952 

0.936 

0.919 

0.904 

0.889 

0.88 

4.  Determination  of  solid  angle:  The  solid  angle,  f)  ,  subtended 
at  the  ground  by  the  bubble  is  porportional  to  the  projected  area  of  the 
bubble  divided  by  Z^  .  Assigning  that  the  bubble  is  spherical,  the  solid 
angle  can  be  represented  by 


n  oc 


Z2 


since  vol  oc  l/p  , 


Q  can  be  expressed  as 
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n  « 


(B-20) 


The  ratio,  O/CIq  ,  of  the  solid  angle  subtended  at  the  ground  by 
the  bubble  at  height  Z  to  the  angle  subtended  by  the  bubble  at  100  meters 
is  tabulated  in  Table  B-IV  for  various  heights  and  lapse  rates,  a  . 


TAfiI£  B-IV 

RATIOS  OF  SOLID  ANGIES 

Z(meters) 


a(“C/lOO  meters) 

100 

200  400  600 

800 

1000 

1200 

1400 

1500 

0 

1 

0.252  0.063  0.029 

0.016 

0.011 

0.007 

0.006 

0.005 

1 

l] 

1 

>  essentially  same  as 

above 

2 

1 

5.  Determination  of  bubble  height  dependence  on  time;  Itois  is  ap¬ 
proached  as  a  buoyant  force  calculation  >diere 

g  =  veight  of  air  , 

oL  SI 

g  =  weight  of  bubble  ,  and 
^a  =  \  ^  • 


Therefore 

(Pa-f’b^^S  =  net  force  on  bubble,  and 

2r 

(Pa-Pb)e^  =  Pb  ^  ^  > 

dt^ 


or 
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(B-21) 


d^Z 

dt^ 


For  an  Isothermal  atmosphere,  substitution  of  (B-15)  and  (B-18)  into 
(B-21)  yields 


(B-22) 


For  finite  lapse  rate,  substitution  of  (B-16)  and  (B-19)  into  (B-21) 

yields 


(B-23) 


The  time  (in  seconds)  required  for  the  bubble  to  rise  from  the  ground 
to  height  Z  is  tabulated  in  Table  B-V  for  various  combinations  of  density 
ratio,  Pt,o/Pao  >  lapse  rate,  a  .  (The  =  0.75  corresponds 

roughly  to  an  initial  temperature  difference  of  100* C.) 


TABLE  B-V 
BUBBIE  RISE  TIMES 


a(*C/lOO  meters) 

^bo'^ao 

0 

200 

0 

0.75 

0 

11.28 

0 

0.65 

0 

15.68 

0 

0.95 

0 

30.1 

1 

0.75 

0 

11.1 

1 

0.85 

0 

15.28 

1 

0.95 

0 

28.35 

2 

0.75 

0 

10.91 

2 

0.85 

0 

14.87 

2 

0.95 

0 

26.8 

Z(meters) 


400 

600 

800 

1000 

1200 

1400 

15.95 

19.5 

22.6 

25.3 

27.8 

30.1 

22.15 

27.00 

31.35 

35.15 

38.6 

41.9 

42.6 

52.4 

60.9 

68.8 

76.8 

83.8 

15.7 

19.2 

22.2 

24.8 

27.15 

29.3 

21.6 

26.5 

30.6 

34.2 

37.4 

40.4 

40.0 

49.0 

56.6 

63.3 

69.4 

75.0 

15.45 

18.89 

21.79 

24.34 

2G.64 

28.7 

21.0 

25.7 

29.65 

33.05 

36.15 

28.95 

37.85 

46.3 

53.1 

58.8 

63.9 

68.5 
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II.  THE  INFRARED  DETECTION  SYSTEM 


A.  Proposed  System 

In  order  to  detect  the  radiation  from  the  gas  in  the  bubble  by  a 
pass!'  e  infrared  search  system,  it  is  necessary  that  the  gas  bubble  emit 
str(^ly  in  the  Infreired  spectrum.  A  number  of  gases  emit  radiation  in  the 
infrared  spectrum,  some  of  \diich  are  constituents  of  the  atmosphere.  A  good 
choice  for  the  heated  gas  of  the  bubble  would  be  a  mixture  of  carbon  dioxide 
and  dry  air.  The  essential  properties  of  the  mixture  eire  that  it  contains  as 
much  carbon  dioxide  as  possible  and  it  must  always  be  less  dense  than  the  at¬ 
mosphere.  The  carbon  dioxide,  and  therefore  the  bubble,  will  radiate  as  a 
black  body  in  narrow  wavelength  intervals  corresponding  to  the  absorption 
wavelengths.  However,  since  it  is  richer  in  carbon  dioxide  than  the  atmos¬ 
phere  and  has  a  higher  partial  pressure  of  CO^,  the  bubble  will  radiate  over 
a  wider  band  of  wavelengths  than  the  atmosphere  due  to  pressure  broadening. 
Ihus,  if  the  detection  system  is  adjusted  to  accept  a  narrow  band  of  wave¬ 
lengths  located  Just  to  one  side  of  the  CO2  absorption  band  of  the  atmosphere 
as  shown  in  Fig.  B-1,  the  CO2  in  the  atmosphere  will  not  absorb  much  of  the 
energy  radiated  by  the  bubble  in  that  narrow  band.  The  atmosphere  will  be 
transparent  to  the  desired  radiation  providing  that  none  of  the  other  constit¬ 
uents  of  the  atmosphew  absorb  energy  in  the  band  of  wavelengths  accepted  by' 
the  detection  system. 

The  choice  of  wavelength  accepted  by  the  detection  system  is  deter¬ 
mined  by  the  absorption  band  of  CO2  which  is  used.  A  wavelength  for  'vdiich  the 
atmosphere  is  transparent  but  is  next  to  a  carbon  dioxide  absorption  wave¬ 
length  would  be  the  most  desirable  wavelength  at  which  to  operate.  A  suitable 
choice  would  be  4.1  microns  which  is  on  the  sharp  edge  of  the  4,3  micron  ab¬ 
sorption  band  of  carbon  dioxide.  The  atmosphere  has  no  constituents  which 
absorb  energy  at  this  wavelength,  consequently  it  is  transparent  at  this 
wavelength. 


Since  4.1  microns  is  a  favorable  wavelength  for  propagation,  the 
next  consideration  is  the  availability  of  suitable  detectors.  In  the  selec¬ 
tion  of  a  suitable  detector,  the  major  factors  to  be  considered  are  sensitivity 
and  optimum  signal -to-noise  ratio.  In  these  considerations,  all  the  incoming 
radiation  is  treated  as  signal  and  the  noise  referred  to  is  that  generated  in 
the  photoconductor.  Available  detectors  have  a  D*  of  about  10^®.  An  estimate 
of  the  signal-to-noise  ratio  can  be  made  by  considering  a.  detector  with  dimen¬ 
sions  1  mm.  x  1  nun.  using  a  90  cycle  per  second  chopper  to  modulate  the  incom¬ 
ing  radiation,  a  5  cycle  amplifier  bandwidth,  smd  liquid  nitrogen  cooling.  A 
signal-to-noise  ratio  of  approximately  55  db.  would  be  obtained  by  this  system. 
This  would  indicate  that  detector  noise  will  be  no  problem  in  this  application. 
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-  6<J 


;,«ospberic  Absorption 


The  significant  noise  in  this  system  will  be  due  to  the  radiation 
from  the  background.  The  system  should  be  sufficiently  flexible  so  that  the 
svin  will  never  shine  directly  on  the  detector  system,  thus  the  noise  radiance 
will  be  from  other  sources.  An  estimate  of  this  noise  Is  0.1  to  2  w/meter^- 
steradlan-mlcron.  The  value  could  be  easily  determined  experimentally. 


B.  Radiance  at  Earth’s  Surface  from  Bubble 


The  Intensity  of  the  radiation  from  the  bubble  of  gas  at  a  tempera¬ 
ture  T  is  given  by  the  equation 


l(\,T)  a  -i  |-  ^  watts/meter^/micron  ,  (B-24) 

x5  [e+Ca/XT.iJ 

where  =  3.74  x  10®  watts  micronVn'Ster^  , 

Cg  =  1.439  X  10^  micron  *K  ,  and 
\  =  wavelength  in  microns. 

The  total  po\«r  emitted  by  the  bubble  is 

W(X,T)  =  l(;>,T)4nr^2  ^  (g.25) 

where  r^jj  is  the  radius  of  the  bubble. 

The  reidlamt  emittance  at  the  earth's  surface  from  a  bubble  of  temper¬ 
ature  T  located  at  a  height  Z  above  the  earth's  surface  is 

l(X,T,r,Z)  =  I(\,T)  r^^/Z^  .  (B-26) 


The  radius  of  the  bubble  is  not  independent  of  the  bubble  temperature 
so  the  r^'-  dependence  in  (B-26)  must  be  removed.  Using  the  perfect  gas  law 
and  the  approximation  h  R  , 
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,  Wb » *bi 
“  L-STiTJ 


(B-27) 


2/3 


SubstitutinG  (B-24)  and  (B-27)  in  (B-2G),  for  4.1  microns,  l(X,T,r,Z)  becomes 


I(X,T,Z) 


3.23  X  10^ 
3510 


2/3 


3  1  JJb 

4  n  pg_ 

L^aJ 


2/3 


1_ 

z2 


(B-20) 


usinc  the  relations  and  Pg  =  ^  Tg  .  Using  on  Mjj  =  30  gm.  , 

p  =  1000  Gr>/m3  gg  average  values,  a  plot  of  the  radiance  at  the  earth's 
surface  from  a  bubble  at  1500  meters  versus  the  difference  between  the  bubble 
temperature  and  the  air  temperature  is  shown  in  Fig.  B-2  for  three  values  of 
the  air  temperature  at  1500  meters.  The  calculations  for  the  data  in  this 
figure  are  given  in  Table  B-VI. 
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TABU)  B-VI 


CALCULATIONS  OF  RALlASiCE  AT  THE  EARTH'S  SURFACE  raCM  A  BUBBU! 
OF  TQffERATURE  Tfe  *  T^+AT  AT  3500  METERS 

ABOVE  THE  EARTH’S  SURF^ 

(X  3  4.1  Microns) 


l(4.1u,T,1500m)  =  5.38x10-3 


1 2/3 


~35To 

.■^-1 


AT  'K 

3510 

'^b 

3510 

."k-i 

/3510 

(TjT^)2/yg  Tt,  _ 

^  l(4.1j.,T^, 1500m) 

For  T. 

®1500m 

=  300* K 

II 

Tg^+AT  (Isothermal) 

10 

11.3 

1.035 

1.023 

80,000 

1.28x10-3 

6.9x10"5 

20 

11.0 

1.07 

1.046 

60,000 

1.78x10-3 

9.6x10"® 

50 

10.0 

1.165 

1.104 

22,000 

5.3x10-5 

2.85x10"'^ 

100 

8.77 

1.33 

1.22 

6,500 

2.05X10"* 

1.1x10”® 

200 

7.0 

1.67 

1.41 

1,100 

1.53x10"^ 

8.2x10"® 

For  T, 

^ISOOm 

-  285*K 

(Adiabatic  Lapse  Rate) 

10 

11.9 

1.035 

1.023 

160,000 

6.4x10“3 

3.44x10"® 

20 

11.5 

1.07 

1.046 

100,000 

1.07x10-3 

5.77x10”® 

50 

10.5 

1.175 

1.114 

37,000 

3.17x10"3 

1.72x10”'^ 

100 

9.1 

1.35 

1.22 

9,000 

1.5x10"* 

8.2x10”^ 

200 

7.25 

1.70 

1.42 

1,400 

1.21x10-3 

6.65x10-5 

For  T„ 

*^1500111 

=  270* K 

(Twice  Adiaoatic 

Lapse  Rate) 

10 

12.5 

1.035 

1.023 

280,000 

3.67x10-5 

1.98x10-5 

20 

12.1 

1.075 

1.05 

180,000 

6x10-5 

3.23x10-® 

50 

11.0 

1.185 

1.12 

60,000 

1.97x10-5 

i.oexio-"^ 

100 

9.5 

1.37 

1.23 

13,000 

1.05x10-* 

5.65x10-7 

200 

7.48 

1.74 

1.45 

1,800 

9.65x10-* 

5.2x10-5 
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III.  EVALUATION  OF  THE  ASSUMPTION  OF  NEGUGIBI£  RADIATION  LOSS 


The  thermod^nasiic  model  used  in  the  preceding  analysis  has  been 
based  on  the  assumption  of  negligible  loss  of  energy  by  thermal  radiation. 

Yet  the  infrared  detection  and  tracking  system  cemnot  operate  without  the  ex¬ 
istence  of  detectable  amounts  of  this  radiation.  The  purpose  of  this  section 
Is  to  examine  this  Inconsistency.  The  approach  tedcen  will  be  to  use  the  re¬ 
sults  of  the  Infrared  detection  analysis  to  determine  the  amount  of  radiation 
required,  and  then  to  determine  whether  or  not  this  amount  of  radiation  would 
be  significant  In  the  thermodynamic  analysis. 

From  the  infrared  analysis,  It  has  been  established  that  a  30-gm. 
bubble,  to  be  detectable  at  a  height  of  1,500  meters,  must  be  some  200*K  above 
the  temperature  of  the  surrounding  air.  It  has  also  been  established  that  the 
bubble  temperature,  in  degrees  Kelvin,  does  not  change  significantly  as  the 
bubble  rises  from  ground  level  up  to  the  1,500-meter  level.  Hence,  a  good 
first  approximation  Is  to  assume  that  the  bubble  temperature  Is  constant  for 
the  calculation  of  energy  loss  by  radiation. 

Assuming  In  addition  that  the  bubble  radiates  as  a  black  body  at  all 
wavelengths,  the  rate  of  loss  of  energy  by  rsidiation  to  the  surroundings  would 
be 


H  =  oAt,(T^-Tj)  ,  (B-29) 

where  =  rate  of  energy  loss,  watts; 

dx 

o  =  the  Stefan-Boltzmemn  constant,  5.67  x  10"®  watts/m^(*K)^; 

A^  =  surface  area  of  the  bubble,  m^; 

T^  =  temperature  of  the  bubble,  “K;  and 

Tg^  -  temperature  of  the  surrounding  air,  *K. 

The  surface  area  of  the  bubble  is  related  to  its  temperature  and  pressure. 

The  radius  of  the  bubble  has  been  computed  from  the  ideal  gas  law  as 
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(B-30) 


»‘b  = 


; 


vhere  >  mass  of  gas  In  the  bubble,  kg.; 

>  universal  gas  constant  7  molecular  weight  of  gas,  13  287 
Joule/kg  *K  if  molecular  weight  =  29;  and 

P-jj  =  pressure  of  gas  in  the  bubble,  newtons  per  m^. 

If  the  bubble  is  assumed  to  be  in  pressure  equilibrium  with  its  surroundings, 
then 


^  ~  ^a  “  ^a^a^a  > 


where 


Pg^  =  air  pressure,  newtons/m^; 

a  air  density,  kg/m^; 

Rg^»  R^;  and 

T  =  air  temperature,  *K. 


Consequently, 


^b  = 


3  MfcTb 


^^Pa^aj 


1/3 


(B-31) 


(B-32) 


8uid  the  surface  area  of  the  bubble  is 


Av  a  4nrv  -  4n 


3  MbTb 


2/3 


(B-33) 
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In  order  to  obtain  numerical  values,  take  the  following  case: 
=  500* K  i 
»  300*K  j 

3  0.030  kg.  >  30  gm.  ;  and 

p„  =  0.05  p_-  (corresponding  to  the  isothermal  case  at 
Z  =  1500  m) 

=  1.09  kg/m3  . 

Then 


dt  " 


(5.67xl0‘®)4n 


3(0.030W 

4n(l.09V( 


500 

300 


)  kg  *K 
)  kg  •k 


2/3 


X  1(500) *-(300)^1  =  1.910  watts  , 

L  -I  m2  Ck)^ 


(B-34) 


The  total  energy  radiated  from  the  bubble  during  its  ascent  is 


^tot  =  11 


(B-35) 


where  t^^  =  ascent  time  to  1,500  meters  (seconds).  Calculations  for  the  rate 
of  rice  of  the  heated  bubble  indicate  that  a  bubble  heated  to  500* K  would  reach 
1,500  meters  in  about  20  sec.  Hence 


^ad  "  1/910  watts  x  20  sec.  =  38,200  Joules  .  (B-36) 


Now  to  ascertain  \diether  this  loss  of  energy  by  radiation  is  sig¬ 
nificant  in  the  thermodynamic  model,  this  energy  may  be  compared  with  the 
total  energy  initially  stored  in  the  bubble.  Assuming  the  bubble  is  heated 
at  constant  volume  to  500®K,  the  total  energy  input  is 
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Q^jj  =  M^CyAT  =  (0.030  kg)  ^716  (200*K)  =  4,300  Joules  (B-37) 


Hence,  it  is  quite  evident  by  ccnparing  and  that  radi¬ 

ation  losses  are  not  negligible.  In  fact,  they  appear  so  large  that  the  sys¬ 
tem  would  be  almost  unworkable  under  the  assumed  conditions.  However,  it  must 
be  pointed  out  that  the  gas  bubble  h\s  been  assumed  to  be  radiating  as  a  black 
body  at  all  wavelengths,  which  is  not  correct  and  results  in  overestimating 
the  energy  loss  by  reidlation.  It  seems  reasonable  to  conclude,  however,  that 
the  heated  bubble  concept  will  not  be  usable  up  to  altitudes  of  1,500  meters. 
If  a  usable  working  range  for  this  system  is  needed,  a  more  detailed  emalysls 
taking  radiation  losses  into  account  in  the  thermodynamic  model  will  be 
required. 
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